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Ifi  spectrum  of  the  Isolated  keto  tautomer  of 
2-axopyrimidioe , its  predicted  spectrum  can  be  compared  vith 
experimental  spectra  of  crystalline  samples  of 
2-oxopyrimldlne  and  matrix  isolated  diketo  tautomer  of 

Simple  statistical  analysis  are  introduced  to 

investigate  the  effect  of  scaling  Che  calculated  force 
fields  in  different  coordinate  systemsr  on  the  predicted  IR 
band  frequencies,  PCDs,  and  intensities  of  these  tautomers. 
The  choice  of  reference  state  in  force  field  and  APTs 

determination  and  the  degree  of  precision  in  ^ initio 
calculations  (convergence  and  integral  threshold)  are 

similarly  studied.  Various  observations  and  recommendations 
are  given.  The  statistical  and  mathemathi cal  approaches  used 
in  studying  these  effects  need  to  be  developed 

substantially . 


The  study  of  pyrimidine  molecules  is  important  in  the 
ongoing  research  on  nucleic  acids.  Nucleic  acid  bases,  which 
are  involved  in  various  key  roles  in  genetic  and  biochenical 
cellular  processes,  are  either  purines  or  pyrimidines.  This 
interpretation  of  vibrational 
.automers  of  Che  2-oxopyrinidine 
l-pyrioidinone  (see  Figure  1-1). 
if  Nonooxopyrimidines 
' theoretical  studies  it  is  worth 


molecule 


the  isolated 
•pyrimidlnol  anc 
1 . 1 Tautonerism 
Before  discussing  tt 


examining  the  tautomeric  behavior  of  Che  system  in  question. 
Nowak  et  al.^  studied  the  IB  and  UV  spectra  of  various 
monooxopyrimidines  in  the  gas  phase,  in  low  temperature 
inert  matrices,  and  in  various  liquid  solvents.  They  found 
that  the  tautomeric  behavior  of  the  2-oxopycimidines  was 
different  from  that  for  the  4-oxopyrimidines . For  the 
4-oxopyrlmidlnes  in  the  gas  phase  and  isolated  in  inert 
matrices  they  found  comparable  concentrations  of  the  keto 
(oxo)  and  enol  (hydroxo)  forms,  whereas  2-oxopyrimidines 
were  found  to  exist  predominantly  in  the  enol  form.  In 
liquid  solvent  systems  both  2-  and  4-oxopyrimidines  existed 
predominantly  in  the  keto  form.  Czecminski  et  al.^  studied 
autoassociaClon  and  CauComerlsm  of  2-oxo-S-x  pyrimidines 
(XaB,Cl,6r)  finding  similar  results  regardless  of  X.  The 


Figure  1-1.  Holeculer  structures  of  2-py rlmidinol , 
2-pyrlsildlnone , cytosine  and  uracil. 


latter  etudiee  were  trade  under  various  condltionsi  in  a 
dilute  Ar  matrix,  more  concentrated  in  an  At  matrix,  in  RBr 
pellets,  as  crystalline  films  (10  and  300°  X),  and  in  an 
amorphous  film  (10°  K).  The  results  from  these  studies  of  IR 
spectra  obtained  for  2'oxopyrlmidine  are  shown  in  Figure  S 
of  Reference  2.  The  results  for  the  2-oxo-5-halogenated 
systems  were  similar.  The  compounds  all  existed  as  enol 
tautomers  when  isolated  in  low  temperature  inert  matrices 
(Figure  5a  of  Reference  2).  This  conclusion  was  evident  from 
the  absence  of  a strong  carbonyl  stretch  absorption  band 
around  1700  cm  and  from  the  appearance  of  strong  bands 
neat  3580  cm"^  where  the  hydroxyl  stretch  is  observed  in 
inert  matrices  and  gaseous  pyrimidines  and  purines . ^ 

When  the  matrix  was  annealed  and/or  enriched,  hydrogen  bond 
association  and  enol  to  keto  tautomertsm  took  place.  This 
conclusion  followed  from  the  appearance  of  C-0  stretch 
absorption  near  1650  cm  ^ as  well  as  from  the  strong  broad 
bands  for  hydrogen  bonded  associates  in  the  3300-2200  cm~^ 
region.  The  RBr  pellets,  crystalline  films  (10  and  300°  K), 
and  amorphous  films  (10°  K1  showed  the  strong  absorption  in 
the  carbonyl  region  neat  1700  cm"^,  indicating  the  keto  form 
as  predominant.  Furberg  and  Solbakk^  determined  from  x-ray 
crystallographic  methods  that  2-oxopyrimldine  existed  in  the 
solid  as  the  keto  form  only,  consistent  with  the  above 
mentioned  results  from  the  IR  spectra,  Czerminski  et  al.^ 
elucidated  a possible  mechanism  for  the  tautomeric 
transition  on  going  from  isolated  molecules  to  form  the 


crystalline  solid. 


In  our  pcesont  study  we  shall  concentrate  on  the 
experinental  results  from  the  dilute  At  matrix  at  10°  R. 
Before  trying  to  understand  more  complicated  matters  such  as 
tautomeric  transitions,  media  effects,  and  perturbations  due 
to  molecular  aseociations,  we  need  full  understanding  of  the 
spectra  of  each  particular  tautomeric  fots  of  the  isolated 
monomers.  The  study  of  2-oxopyrimidine  In  the  dilute  Ar 
matrix  is  such  a case,  since  only  the  enol  form  occurs,  in 
this  sudy  we  shall  work  with  a more  detailed  and  better 
resolved  and  still  unpublished  spectrum  of  2-oxopyrlmidine 
in  the  dilute  Ar  matrix  than  is  shown  in  Figure  5a  of 
Reference  2.  This  spectrum  has  been  provided  to  us  by  the 
Institute  of  Physics  of  the  Polish  Academy  of  Sciences.^ 

We  will  carry  out  theoretical  calculations  of  the  IB 
spectrum  for  each  one  of  the  tautomeric  species  of  the 
2~  oxopyrimidine  molecule.  The  calculations  will  be  compared 
mainly  with  the  above  mentioned  experimental  spectcum  of 
2-oxopyrlmidine  in  an  Ar  matrix. 

1 . 2 Theoretical  Perspective 

Quantum  mechanical  calculations  utilize  extensively 
molecular  orbital  methods  (HO).  The  two  major  approaches  may 
be  classified  as  semiempirical  and  ab  initio  methods.  The 
first  may  utilise  a simplified  Hamiltonian,  but  is 
characterised  by  the  use  of  experimental  results  to  evaluate 
parameters  used  in  the  calculations.  The  second  one  uses  Che 
correct  Hamiltonian  and  evaluates  all  parameters  by 


reatheinathical  methods,  avoiding  use  of  experimental  results. 
We  will  be  concerned  here  with  ab  initio  studies  using 
Hartree-Fock  iHF)  self-consistent  field  fscF)  calculations. 

In  this  approach  the  wavefunctlon  of  the  system  is 
written  as  the  Slater  determinant  formed  from  a set  of 
occupied  spin-orbitals.  These  spin-orbitals  are  just 
products  of  ordinary  molecular  orbitals  inO)  and  spin 
functions  for  each  electron.  The  molecular  orbitals  are 
themselves  written  as  linear  combinations  of  atomic  orbitals 
(AO).  Each  atomic  orbital  is  expressed  as  a linear 
combination  of  simple  basis  functions.  The  energy 


adjusting  coefficients  in  these  linear  combinations  so  as  to 
obtain  a HP  SCF  wavefunction  for  the  electronic  ground  state 
for  a particular  nuclear  arrangement. 

various  basis  functions  have  been  used  in  the  past. 
The  earliest  ones  were  Slater  type  orbitals  (5TO)  used  as 
approximations  to  atonic  orbitals.  The  HO  SCF  calculations 
use  either  minimal  or  extended  basis  sets.  In  theory  an 
infinite  number  of  AOs  can  be  used  to  construct  a mo,  but 
during  calculations  approximations  have  to  be  made.  A 
minimal  basis  set  of  AOs  consists  only  of  occupied  inner  and 
valence  shells,  whereas  an  extended  basis  set  also  includes 
unoccupied  higher  shell  AOs. 

Because  of  computational  difficulties  with  STOs  (some 
integrals  cannot  be  readily  evaluated  by  the  computer),  most 
worlcers  have  replaced  them  with  linear  combinations  of 
Gaussian  functions.^  Among  the  minimal  basis  sets  used  one 


find  the 


review  by  Pople' 


stated  that  results  using  5TO-3G  and  ST0-4G  basis  sets 
agreed  quite  well  with  calculations  based  on  exact  STOs 
studies.  A better  description  of  the  wavefunctlon  is 
provided  by  a linear  conbination  of  Gaussian  functions 
called  split  valence  basis  sets.  The  number  of  functions  is 
doubled  for  the  valence  shell  to  provide  extra  flexibility 
to  the  molecular  orbitals.  Among  these  one  csn  usually  find 
the  3-21c,^^  4-21,^^  and  4-31C^^  basis  sets  as  comnoniy  used 
ones.  The  extended  basis  sets  are  those  including 
polarization  functions  in  the  description  of  the 
wavefunctlon.  Commonly  used  extended  basis  sets^^'^®  are 


increased  the  predictions  from  this  wavefunctlon  will  be 
better.  It  is  expected  that  additional  improvements  will 
result  from  including  corrections  to  the  wavefunctlon  due  to 
configuration  interaction,  i.e,,  linear  combinations  of 
Slater  determinants.  As  far  as  geometry  optimizations  and 
force  field  evaluations  are  concerned,  these  expectations 
are  generally  fulfilled.  But  the  high  cost  and  complexity  of 
these  calculations  on  polyatomic  molecules  have  limited  such 
studies.  It  is  beyond  the  scope  of  the  present  work  to 
evaluate  basis  sets  and  correlation  effects.  The  reader  may 
find  several  reviews  which  deal  with  these  aspects  in  more 
8,17-19 


depth,* 


In  the  present  work  we  shall  look  almost  exclusively 
at  results  froin  studies  with  small  basis  sets  such  as  ST0-3Q 
and  4-31G  at  the  HF  SCF  level  connected  with  the  prediction 
and  interpretation  of  IFt  spectra  of  polyatomic  molecules.  We 
have  classified  these  results  as  geometry  optimization, 
choice  of  reference  state,  determination  of  frequencies, 
scaling  of  force  fields,  and  determination  of  intensities. 
These  subjects  are  discussed  in  more  detail  in  the  following 
sections . 

1.2.1  Geometry  Optimization 

determination  of  the  geometry  of  a chemical  system  at  the 
minimum  of  its  potential  energy  curve.  A given  point  on  this 
potential  curve  or  surface  is  the  sum  of  the  electronic 
energy  and  of  the  nuclear  repulsion  energy  at  a given 
nuclear  arrangement  of  the  system.  In  the  past  the  points 
along  these  curves  were  determined  from  calculations  of  the 
energy  made  at  each  separate  point  and  the  equilibrium 
geometry  was  determined  by  fitting  them  to  given  functions. 
These  explorations  were  carried  out  partially  by  varying 
certain  parameters  while  holding  others  constant,  with 
recent  improvements  in  programming  direct  analytical 
calculations  of  energy  derivatives  have  been  carried  out,^® 
so  that  the  energy  of  the  system  can  be  minimized  with 
respect  to  all  variables  simultaneously.  Thus  quite  valuable 
theoretical  equilibrium  geometries  have  been  determined. 


PQpIe°  reviewed  at  some  length  geometry  predictions  at 
the  Hartree-Pock  level.  He  pointed  out  that  at  the  tine  of 
his  review  nost  work  was  done  at  this  level.  For  large 
raolecules  {around  10  atoms}  most  studies  have  been  done  with 
minimal  basis  sets.  For  smaller  molecules  larger  basis  sets 
have  been  used.  The  effects  of  higher  levels  of 

approximation  such  as  configuration  interaction  have  been 
studied  to  a very  small  extent  for  large  molecules.  It  is 
claimed  Chat  HF  theory  using  minimal  basis  sets  {such  as 
STO-3G)  has  been  moderately  successful  in  determining 

equilibrium  geometries  for  closed  shell  systems.  Bond 
lengths  and  angles  are  usually  given  with  an  accuracy  of 
0.03  A and  respectively.  Use  of  split  valence  basis  sets 
(such  as  4-31G)  improved  predictions  of  bond  length  but  not 
bond  and  dihedral  angles.  Polarization  functions  added  to 
split  valence  basis  sets  (such  as  S-31G  ) improved  the 
prediction  of  the  angles  and  gave  equilibrium  geometries 
close  to  the  HF  limit.  The  bond  lengths  were  systematically 
underestimated  by  0.01-0.02  A.  Very  few  predicted  geometries 
with  corrections  for  electron  correlation  have  been  studied. 
It  seems  that  no  major  improvements  are  obtained  as 
compared  to  the  above  mentioned  basis  sets. 

1.2.2  Choice  of  Reference  State 

The  potential  energy,  V,  of  a chemical  system  can  be 
written  as  a power  series  expansion  in  the  displacement 


coo rdinates 


V(xi,*j,  . . , ,x3K|  - * E(  SV/a*i)g'Axl  t 

0.5*  E(  3^V/3xl3xj  ^j'fixiAxj  ♦ ....  (1-1) 

where  xl,x2, x3K  ace  the  cartesian  coocdinaCes  cf  Che  h 

nuclei.  The  first  terra,  V , is  Che  potential  energy  of  the 
nolecule  at  eguilibriuin.  The  second  terra  vanishes  at  the 

in  the  harraonic  approximation  since  higher  order  derivatives 
are  neglected.  Thus 

V(xi,xj x3H)-0.5*EFi j'flxiflxj  (1-21 

where  Fi)  ( > 3^v/3xi 3xj ) are  known  as  the  force  constants. 
This  appcoxiiiatiDn  for  the  potential  energy  function  is 
acceptable  for  small  displacements  near  Che  equilibriura 
position  as  in  low  amplitude  vibrations. 

From  the  force  constants  the  IP  frequencies  can  be 
obtained.  In  theoretical  calculations  a given  molecule  will 
be  distorted  slightly  frora  the  equilibriura  configuration  in 
order  to  determine  the  force  constants.  This  implies  a 
priori  knowledge  of  the  eguilibcium  configuration.  The 
choice  of  this  configuration  that  we  can  call  reference 
state  is  a controversial  topic.  Different  criteria  ate  used 
to  select  it.  Several  typical  examples  are  cited  below. 


Some  workecs  have  utilized  experimental  equilibrium 
geometries  as  ceCerence  states . In  some  cases  better 
force  constants  were  obtained  than  with  theoretically 
determined  reference  states.  Usage  of  experimental 
geometries  shortened  the  computations  since  there  was  no 
need  for  geometry  optimization  calculations. 

Others  preferred  theoretical  reference  states  and  the 
choices  were  various.  Pulay  et  al.,^^  and  Blom  and  Altona^*^ 
have  used  empirically  corrected  SCF  geometries.  They 
evaluated  the  equilibrium  geometry  with  a given  basis  set 
and  then  utilised  empirical  corrections  to  approximate  it  to 
experimental  observations.  They  claimed  that  the 
experimental  geometry  is  suspected  of  giving  random  errors 
in  the  evaluation  of  force  constants,  experimental 
geometries  could  be  different  depending  on  the  experimental 
technique  employed.  Further  still  most  of  them  are 
vibrational  averages  which  are  different  from  true 
equilibrium  geometries.  This  randomness  is  reflected 
significantly  in  the  evaluated  force  constants.  Fogarasi  and 
Pulay^^  pointed  out  that  for  email  molecules  it  is  better  to 
use  high  quality  theoretical  calculations  for  the 
determination  of  the  equilibrium  geometry.  Then  determine 
the  force  constants 

at  a lower  theoretical  level,  e.g.,  SCF,  as  demonstrated  by 
the  work  of  Meyer  and  Rosmus  on  diatonic  hydrides. 

Another  group  of  researchers  used  theoretically 
deternined  equilibrium  geometries  as  reference  states 
without  any  modifications  at  all.^^'^^  They  evaluated  the 


theoretical  ground  state  with  a given  level  of  theory,  e.g., 
a given  basis  set.  Then  they  determined  the  force  field  at 
the  same  level  of  theory  using  this  geometry.  Thus  we  have 
examples  such  as  the  worh  by  Nishiraura  et  al.^^  on  uracil. 
They  optimised  uracil  at  the  HF  SCF  level  with  an  ST0-3G 
basis  set.  Then  they  determined  Che  force  field  at  the  same 
level  of  theory  with  this  geometry.  Pople  went  as  far  as 
saying  that  it  mates  no  sense  to  depart  from  this 

approach . These  workers  claim  that  in  this  way  a more 

systematic  approach  of  vibrational  studies  can  be  carried 
out,  and  data  could  be  easily  compared  between  different 
groups.  This  is  akin  to  Che  concept  of  "theoretical  model 

chemistry"  which  uses  a uniform  level  of  theory 

throughout . In  a recent  review  by  Hess  and  Schaad^^  on 
ab  ini  tip  calculations  of  vibrational  spectra  they  pointed 
out  that  this  procedure  is  at  least  equivalent  to  the  other 
mentioned  choices  of  reference  state. 

1.2.3  Determination  of  Vibrational  Frequencies 

The  vibrational  frequencies  are  determined  from  Che 
force  constants  (Equation  1-2)  at  the  harmonic 

approximation.  The  force  constants  are  determined  in  turn  by 
numerical  or  analytical  gradients  of  the  potential 

energy. Most  of  Che  work  on  medium  to  large  polyatomic 
systems  has  been  carried  out  at  the  Hartree-Fock  SCF  level 
with  minimal  basis  sets.  This  is 
cost  of  computer  time. 


mainly 


These  calculations  usually  overestimate  the  force 
constants  and  frequencies  when  compared  to  experimental 
results.  Hess  and  Schaad,^^  reviewing  an  impressively  large 
number  of  molecules,  concluded  that  frequency  errors  are 
usually  in  the  5-lDt  range  for  Hartree-Pock  SCP  calculations 
with  minimal  basis  sets.  Pople  et  al.^^  carried  out  HP  SCP 
calculations  with  3-21G  and  6-31G  basis  sets  on  37  email  to 
medium  sieed  molecules  13-10  first  row  atoms).  They 
concluded  that  theoretical  frequencies  exceeded  experimental 
ones  by  ISO  cm'^  (about  121)  at  the  3-21G  level.  Expansion 
of  the  basis  set  to  6-31G  (including  d functions  on  atoms 
other  than  hydrogen)  did  not  improve  agreement  with 
experiment.  Other  researchers  concluded  that  usage  of  double 
seta  or  split  valence  basis  sets  give  frequencies  that  are 
larger  than  experimental  ones  by  j.jsv . ■ 23 ■ 29 , 35-37 

Pogarasi  and  Pulay  provided  us  with  detailed 
observations  on  the  results  from  calculations  with  minimal 
basis  sets  (5T0-36)  and  medium  sized  ones  (3-21G,  4-31G, 
4-21)  for  force  constants  of  medium  to  large  molecules.  They 
pointed  out  that  diagonal  terms  of  force  fields  from  minimal 
basis  set  calculations  were  much  overestimated  when  compared 
with  medium  sized  basis  sets.  The  STO-3G  off-diagonal  force 
constants  showed  random  variations  whereas  medium  sized 
basis  sets  cases  were  more  systematic.  For  medium  sized 
basis  sets  they  found  that  diagonal  stretching  and  bending 
force  constants  were  systematically  overestimated  by  10-154 
and  20-30%  respectively.  For  off-diagonal  force  constants 


the  deviations  were  lees  systenatici  large  force  constants 
were  in  error  by  10-301 « and  small  ones  had  absolute  errors 
in  the  0.05-0.10  aJ/A"  values. 

1.2.4  Scaling  of  force  Fields 

As  seen  in  the  previous  section  H?  SCP  calculations  on 
medium  to  large  molecules  give  frequencies  that  are  usually 
overestimated  when  compared  to  experimental  values.  This 
discrepancy  is  usually  attributed  to  the  usage  of  snail 
basis  sets  and  lack  of  electron  correlation.  This  attitude 
is  in  turn  due  to  computer  limitations  and  high  cost. 
Because  it  is  not  possible  to  do  exact  calculations  on 
medium  to  large  chemical  molecules  many  researchers  have 
systematically  determined  force  fields  at  the  HP  SCP  level 
with  small  basis  sets,  and  then  scaled  them  eo  as  to 
approximate  them  to  the  "true  experimental  force  field." 
various  scaling  approaches  have  been  carried  out. 

A theoretical  force  field  can  be  visualized  as  a 
matrix  containing  diagonal  and  off-diagonal  values  of  force 
constants  in  a given  coordinate  system.  The  most  usual 
coordinate  system  in  scaling  procedures  is  the  internal 
coordinate  set.  e.g.,  bond  stretches,  angular  deformations, 
etc.  The  researcher  can  vary  any  one  of  the  matrix  elements 
so  as  to  fit  his  theoretical  force  field  to  experimental 
observations.  Two  major  types  of  variations  have  been 
carried  out:  diagonal  elements  only,  and  all  elements  of  the 


Among 


examples  in  which  variation  only  of  diagonal 


elements  were  considered  we  find  several  interesting 
examples.  Pulay  and  Meyer^*^  scaled  the  theoretical  force 
fields  of  ethane,  ethylene,  and  acetylene.  They  reduced  all 
diagonal  stretching  force  constants  by  lOt,  and  those  for 
all  bendings  by  20t.  With  this  simple  adjustment  they 
reproduced  experimental  frequencies  within  20-30  cm 
{1. 0-5.9%  accuracy).  Other  researchers^^  have  fitted  the 
diagonal  elements  of  the  force  field  to  experimental 
frequencies  with  excellent  results.  Fogarasi  and  Pulay^^ 
pointed  out  several  disadvantages  with  these  scaling 
procedures.  They  claimed  that  unless  the  calculated 
off-diagonal  force  constants  are  grossly  in  error,  different 
force  fields  can  be  made  to  fit  experimental  frequencies 
analogously.  They  said  that  in  the  study  of  formaldoxime  by 
Pouchan  et  al.,^^  the  same  agreement  was  obtained  with 
experiment  by  scaling  STO-3G  and  S-31G  force  fields  in  spite 
of  these  two  being  significantly  different.  In  this  way  the 
relative  values  of  diagonal  theoretical  force  constants  were 
lost  completely,  and  mi ssassignraents  of  experimental  bands 
could  not  be  detected. 

Scaling  of  all  the  elements  of  the  force  field  seems  a 
more  logical  approach.  Blom  and  Altona^*^  scaled  the  4-31G 
force  fields  of  ethane,  propane,  and  cyclopropane.  They 
separated  the  diagonal  force  constants  in  the  following 
groups:  C-C  stretching,  C-H  stretching,  bendings  and 

roctcings,  torsions,  CH2  wagging  (in  cyclopropane).  All  the 
were  placed  in  one  group.  Then  by 


f-diagonal  elements 


they  found  scaling 


standard  Iterative  procedures'^ 
factors  to  oiulCiply  these  groups  by  so  as  to  obtain  best  fit 
with  experimental  results,  zn  this  way  they  came  up  with 

constants.  They  were  able  to  reproduce  experimental 

freguBocies  to  about  10  cm'^  (1%  accucacyl.  various  other 
calculations  along  these  lines  gave  similar  results. 

Further  scaling  can  be  carried  out  by  finding  more 
conversion  constants.  The  work  by  wiberg  et  al.^^  on 
pyridine  corrected  all  the  4-31G  force  field  diagonal 
elements  with  individual  scaling  factors.  The  off-diagonal 
elements  were  divided  into  four  groups  {one  for  each  of  the 
symmetry  groups  of  pyridine),  and  a single  scaling 

constant  for  each  different  group.  The  following  equation 
{Equation  1-3)  shows  how  Pulay  and  cowor)tBrs^^ ' scaled 
their  theoretical  force  field,  by 


(1-3) 


where  the  conversion  constants  of  the  diagonal  matrix  C gave 
the  force  field  F which  best  fitted  experimental  results.  Zn 
this  way  the  diagonal  elements  of  F^^  such  as  Fii  and  Fjj 
were  multiplied  by  Cii  and  CJj  respectively.  Off-diagonal 
elements  such  as  Fi}  were  multiplied  by  the  geometric  mean, 

( Ci i*C  j j ) *^  ’ ^ . fogarasi  and  Pulay^^  argued  that  this  approach 
was  more  physically  significant  than  conversions  of  the  type 
made  by  Blora  and  Altona^^  which  made  force  constants 
coupling  different  coordinates  share  the  same  scaling 
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constants.  Hlpps  and  Poshusta^*^  have  utilised  the  arithmetic 
mean  instead  of  the  georaetric  one  discussed  here  for  scaling 
coupling  force  constants.  A good  reproduction  of 
experimental  frequencies  has  been  obtained  by  these 
approaches . 

As  we  have  seen  conversion  constants  have  been 
utilised  to  obtain  the  best  fit  between  theoretical  and 
experimental  results.  A different  philosophy  of  scaling  has 
been  adopted  by  other  investigators^^ ' which  we  will 
shall  call  upgrading.  In  this  case  conversion  constants  are 
utilized  to  convert  a lower  level  theoretical  force  field 
into  a higher  level  one.  for  example  Nishimura  et  al.^*^ 
determined  the  force  field  of  uracil  at  the  HP  8CF  level 
with  an  STO-3G  basis  set.  Then  they  utilized  4-31G/STO-3G 
conversion  constants  to  upgrade  the  ST0-3G  force  field  into 
one  which  should  reaerahle  a force  field  evaluated  at  the  HP 
SCF  level  with  a 4-31G  basis  set.  The  conversion  was  done  by 
Equation  1-3  where  was  the  STO-3G  force  field,  and  F was 
the  upgraded  or  "near  4-31G"  one.  The  conversion  constants  C 
were  obtained  from  comparisons  of  previous  calculations  at 
the  5TO-3G  and  4-31G  levels  for  urea,  acrolein,  and 
formamide,^®’®®  It  was  emphasized  that  no  attempts  were  made 
to  fit  theoretical  results  to  experimental  observations. 
There  are  various  reasons  for  applying  this  upgrading 
technique,  primarily  the  cost  of  the  calculations.  As  the 
size  of  the  molecule  increases  so  does  the  complexity  of  the 
calculation  and  the  cost.  By  upgrading  low  level  force 
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fields  such  as  those  from  ST0-3G  calculations,  better 
quality  force  fields  can  be  obtained  at  reasonable  prices. 
Another  inportanc  reason  to  use  this  approach  Is  to  know  how 
far  ^ initio  calculations  can  help  in  understanding 
chemical  systems  without  referring  to  experimental 
observations.  In  this  fashion  purely  theoretical  predictions 
may  be  utilised  to  deduce  experimental  results. 

1.2.5  Peternination  of  vibrational  intensities 

Predicted  IR  spectra  should  include  the  intensities  of 
the  bands  as  well  as  the  frequencies  and  normal  modes, 
unluckily  there  is  less  knowledge  available  about 
intensities  than  on  force  constants.  This  is  quite  evident 
on  studies  at  the  HF  5CF  level  on  medium  to  large  polyatomic 
systems . Fogarasi  and  Pulay^®  stated  that,  from  the 
data  available,  fairly  large  basis  sets.  Including 
polarisation  functions,  were  necessary  to  obtain  reasonable 
values  for  predicted  Intensities  (reliable  to  within  50%). 
They  further  stated  that  calculations  using  smaller  basis 
sets  can  be  considered  to  give  correct  qualitative 
predictions;  i.e.,  they  could  differentiate  between  weak, 
medium,  and  strongly  absorbing  bands. 

In  actual  calculations  of  band  intensities,  the 
derivatives  of  the  dipole  moment  of  a molecule  with  respect 
to  a coordinate  system  are  evaluated  first.  These  dipole 
moment  derivatives  are  then  transformed  into  derivatives 
with  respect  to  normal  modes  via  a transformation  matrix 
(the  L matrix,  see  Chapter  2).  This 


transformation  matrix 
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turn  is  obtained  fcoo  the  dlagonallzation  of  the  calculated 
force  field.  Thus  it  can  be  said  that  the  values  of  the 
calculated  band  intensities  depend  on  the  interaction 
between  Che  calculated  dipole  nonent  derivatives  and  the 
calculated  force  field.  Person  reviewed  these  tranfornations 
in  great  detail. 

Hehre  et  al.^^  compared  ab  initio  calculations  of 
electric  dipole  moments  with  especimental  ones.  They  focused 
mainly  on  comparisons  of  magnitudes  since  experimental  data 
on  signs  and  directions  of  dipole  moments  are  scarce  and 
have  considerable  uncertainty.  They  considered  studies  on  a 
series  of  diatonic,  small  polyatomic,  hypervalent, 
hydrocarbon,  and  monosubstituted  benzene  molecules  at  the  HF 
SCF  level  with  ST0-3G,  3-21G,  3-21G*,  and  $-31G*  basis  sets. 
In  general  they  claimed  moderate  success  with  3-21G  and 
6-31G  basis  sets,  and  recommended  caution  with  5TO-3G  basis 
set  results.  Among  their  results  they  pointed  out  that  for 
diatomic  and  small  polyatomic  systems  the  STO-3C,  3-21G, 
3-21G  , and  6-31G  basis  sets  gave  calculated  dipole  moments 
larger  than  experimental  ones  by  0.9,  0,49,  0.34,  and  0.30  b 
respectively.  For  larger  compounds  containing  heteroatoms 
these  errors  in  calculated  dipole  moments  for  STO-3G,  3-21G, 

0.30  0 respectively.  For  hydrocarbon  compounds  these 

deviations  for  STO-3G,  3-21G,  and  6-31G  basis  sets  were 
0.10,  0.06,  and  0.07  D respectively.  For  monosubstituted 
benzenes  the  deviation  for  5T0-3G  basis  set  was  0.5  I>. 


The  effect  of  electron  correlation  on  Che  calculation 
of  dipole  moirients  has  been  briefly  considered.  As 
discussed  above  it  vas  indicaCed^^  Chat  near  Che  HF  limit 
Che  magnitudes  of  electric  dipole  momente  were  in  error  by  a 
few  tenths  of  a debye.  This  discrepancy  was  significant 
since  magnitudes  of  most  dipole  moments  are  within  0-f  D. 
This  error  then  was  particularly  significant  for  molecules 
with  very  small  dipole  moments  such  as  carbon  monoxide.  The 
experimental  dipole  moment  of  CO  is  0.11  D in  the  direction 
■ceo*  . The  dipole  moment  calculated  at  near  HP  limit  was 
0.21  D in  the  '*'CeO  direction.  Inclusion  of  electron 
correlation  in  the  calculations  corrected  the  near  UP 
prediction  of  the  magnitude  and  predicted  the  sign  of  the 
dipole  moment  in  CO  in  the  right  direction. 

It  is  clear  that  the  above  mentioned  observations  on 
calculated  dipole  moments  provide  us  only  with  approximate 
trends,  note  work  is  nesded  in  order  to  have  more  precise 
knowledge.  The  calculation  of  dipole  moments  is  subject  also 
to  a decision  on  the  reference  state.  This  has  been  pointed 
out  by  uehre  et  al.^^  who  used  experimental  or  calculated 
geometries  on  different  occasions.  They  pointed  out  that  in 
order  to  evaluate  the  reproduction  of  experimental  dipole 
moments  by  different  basis  sets,  two  sources  of  error  should 
be  distinguihedr  limitations  of  the  theory  in  describing 
charge  distribution,  and  improper  choice  of  molecular 
geometry . 


calculations 


intensities  is  concerned,  the  situation  is  more  complex  than 
for  dipole  moment  cases.  As  previously  pointed  out,  little 
Is  known  about  ^ initio  calculations  of  tR  intensities  Cor 
large  molecules.  From  the  conclusions  on  dipole  moments  it 
becomes  clear  that  future  studies  need  to  consider  how  the 
differences  between  calculated  dipole  moment  derivatives  and 
experimental  ones  differ  from  corresponding  differences  for 
dipole  moments.  Further  still  Che  choice  of  reference  state 
and  the  interplay  between  these  derivatives  and  force  fields 
has  Co  be  faced  in  the  light  of  band  intensities 
observations . 


The  primary  goal  of  this  study  is  to  interpret  the 
experimental  IR  spectrum  of  2-oxopyriraidine  measured  in  a 
dilute  Ar  matrix  at  10^  K,  In  order  to  do  so  reliable 
spectra  for  2-pyrimidinol  and  2-pyrimidinone  have  been 
calculated.  The  decision  of  the  level  of  theory  to  be  used 
in  these  calculations  was  based  on  the  limitations  of 
available  knowledge  and  programming  costs.  Secondary  goals 
are  considerations  about  Che  effect  of  different  coordinate 
systems,  convergence,  integral  threshold,  and  choica  of 
reference  state  in  Che  acaling  {upgrading)  of  force  fields, 
and  the  nature  of  predicted  IR  spectra  I frequencies,  normal 
modes,  and  intensities). 

Chapter  2 summarises  Che  theory  used  in  ab  initio 
calculations,  and  the  usage  of  the  latter  In  the 
determination  of  frequencies,  normal  modes,  and  intensities 


o£  IR  fundamental  bands  for  the  present  study.  Chapter  3 is 
concerned  with  an  accurate  calculation  of  the  IR  spectrum  of 
2-pyrinildlnol . This  spectrum  Is  used  to  aesign  the 
experimental  fundamental  bands  of  2-pyrimidinol  in  the  low 
temperature  dilute  Ar  matrix  IR  spectrum  of  2-oxopyrimidine . 

Chapters  4 and  S are  concerned  with  the  effect  of 
coordinate  system,  convergence,  integral  threshold,  and 
choice  of  reference  state  on  the  scaling  of  force  fields  and 
Che  nature  of  calculated  IR  frequencies,  normal  modes  and 
intensities.  Simple  statistical  approaches  have  been 
introduced  to  evaluate  the  quality  of  the  calculated  IR 
spectrum.  Chapter  6 is  concerned  with  the  calculation  of  the 
IR  spectrum  of  2-pytiraidinone . Some  calculations  for  uracil 
have  been  considered  in  order  to  help  in  these 
considerations. 

As  a more  generaliaed  aim  the  results  of  the  present 
work  will  be  used  in  the  future  for  the  interpretation  of  IR 
spectra  of  halogenated  oxopyrimldines  and  of  spectra  of 
nucleic  bases,  particularly  cytosine  (figure  1-1).  The 
advantage  of  examining  the  IR  spectrum  of  the  low 
temperature  dilute  Ar  matrix  of  the  2-oxopyrlmidine  molecule 
is  that  it  provides  a rare  opportunity  of  studying 
exclusively  the  IR  spectrum  of  the  hydtoxo  (enoU  tautomer 
of  2-oxopyrimidine  (i.e.  2-py timldinol ) . This  in  turn  can  be 
utilised  in  helping  with  the  interpretation  of  the  IR 
spectra  of  nucleic  bases.  Spectra  of  these  compounds  are 
often  made  up  of  contributions  from  different  tautomeric 
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forms  of  the  nucleic  base.  A case  in  point  is  the  spectrum 
of  cytosine.  One  of  the  tautomeric  forms  of  this  nucleic 
base  is  very  similar  to  the  structure  of  2-pyrimidinol . 

The  statistical  approaches  introduced  in  Chapters  A 
and  5 are  given  with  Che  intent  of  emphasising  that 
evaluation  of  the  quality  of  a predicted  IB  spectrum  in 
terms  only  of  band  frequencies  is  not  necessarily 
sufficient.  The  observations  are  understood  as  simple  trends 
and  should  be  taken  as  qualitative  only.  In  future  work 
there  should  be  more  concern  in  considering  simultaneously 
the  predicted  frequencies,  normal  modes,  and  Intensities  of 


CHAPTER  2 
THEOBY 

In  this  chapter  we  shall  outline  the  nethode  of 
calculation  used  for  the  prediction  and  Interpretation  of  IB 
spectra.  He  shall  divide  the  discussion  into  two  parts:  ab 
initio  calculations  and  vibrational  spectra.  The  first  one 
considers  the  ^ initio  molecular  orbital  calculations,  and 
Che  second  one,  Che  usage  of  such  results  in  Che  calculation 
of  Che  frequencies,  normal  modes,  and  intensities  of  the 
vibrational  spectra. 

2.1  Ab  Initio  Calculations 

All  the  calculations  were  done  by  the  restricted 
Hartree-Foch  self-consistent  field  raethcd  (RHf  SCF)  with 
STO-3G  and  4-31G  basis  sets.  The  molecules  under 
consideration  were  2-pyrimidinol  and  2-pyriDidlnone . The 
programs  used  were  GAUSSIAN  and  GAUSSIAN  82.^^  GAUSSIAN 

76^^  was  used  with  Che  Amdahl  470  computer  at  the  North  East 
Regional  Data  Center  (NEBDCI,  and  GAUSSIAN  82^^  with  the 
VAX-780  of  the  Quantum  Theory  Project  (QTP),  Department  of 
Chemistry  (University  of  Floridal.  These  programs  were  used 
at  various  times  for  different  purposes.  The  reasons  for  the 
choice  of  one  over  the  other  were  due  to  factors  such  as 
cost,  time  of  installation,  relation  to  time  of  computation, 
and  nature  of  Che  computations.  This  section  is  organized 
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in  two  pacts:  GAUSSIAN  calculations  and  GAUSSIAN  76^^ 

calculations . 

2.1.1  GAUSSIAN  82  Calculations 

Various  types  of  calculations  were  done  with  the 
GAUSSIAN  82  pcogcam.^^  They  ace  divided  under  three 
headings:  geometry  optimizations,  RHF  5CF  4-31G 

calculations,  and  SHF  SCF  STO-3G  calculations.  The  first  one 
is  concerned  with  the  determination  of  the  equilibrium 
geometries  of  both  tautomers  of  2-oxopyrimidine . The  second 
one  loots  at  the  determination  of  the  force  field  and  atomic 
polar  tensors  of  2-pyr imidinol  with  a 4-31G  basis  set.  The 
third  one  considers  the  determination  of  the  force  fields  of 

both  tautomers  of  2-oxopyrimidine  with  an  STO-3C  basis  set. 

2. 1.1.1  Geometry  optimisations 

The  SFGS  optimization  procedure  described  by  Head  and 
Zerner^^  (installed  as  a GAUSSIAN  82^^  subroutine)  was  used 
with  its  default  parameters  for  all  optimizations.  This 
subroutine  performs  simultaneous  optimisation  of  all 
variables  by  analytic  gradient  methods.  Ail  calculations 
were  done  at  the  RHF  SCF  level. 


2-oxopyrimidine  were  determined  with  an  STO-3G  basis  set. 
The  equilibrium  geometry  of  2-pycimidinol  was  also 
determined  with  a 4-31G  basis  set.  The  Initial  guesses  for 
the  ST0-3G  optimization  of  both  tautomers  were  those 
published  from  the  MINDO/3  calculations  by  Czerminski  et 
al.^^  The  STO-3G  optimized  geometries  are  given  in  Chapters 
4 and  8.  The  initial  guess  for  the  4-31G  optimization  of 


2-pycioidinol  wa£  the  geometry  optimized  with  the  STO-3G 
basis  set.  The  4-31G  optimized  geometry  of  2-pycioidinol  is 
given  in  chapter  3. 

2. 1.1. 2 RHF  5CF  4-31G  calculations 

The  subroutine  FORCE  of  GAUSSIAN  82^^  was  used  to 
determine  the  cartesian  force  field  and  the  atomic  polar 
tensors  Iaft)  of  2-pyrimidinoX.  The  calculation  was  done  at 

parameters  of  GAUSSIAN  82^^  were  used  (including  convergence 
• 10'®  and  integral  threshold  ■ 10'^®). 

This  subroutine  was  done  a total  of  86  tines,  one  was 
for  the  reference  state  which  was  at  the  above  mentioned 
4-31G  optimized  geometry  (Section  2. 1.1.1),  and  the 
remaining  58  were  done  for  each  of  the  displaced  geometries. 
Each  atom  of  the  reference  state  was  displaced  one  at  a time 
in  the  +x.-x,+y,-y,  and  *z  direction  by  0.005  A.  The 
cartesian  forces  and  the  dipole  moments  were  calculated  for 
the  reference  state  and  each  of  the  displaced  geometries. 

Each  run  of  FORCE  subroutine  yielded  a 1 x 33  matrix 
of  cartesian  forces  acting  on  each  atom.  An  element  of  this 
matrix  can  be  called  f^^  where  i stands  for  a given  atom, 
j for  one  of  the  three  cartesian  coordinates  (x,y,z),  and  t 
for  a positive  displacement.  Similarly  f^j  stands  for  a -j 
displacement  of  atom  i.  From  these  matrices  a 33  x 33 
cartesian  force  field,  F , was  constructed  by 
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^ij  ■ * 15. 5691/(2(0. 005/0. 529167)) 

(2-1) 

where  wae  an  element  of  matria  In  ndyn/A.  Similarly 
from  the  dipole  moment  conponenta  of  each  run  of  FORCE 
subroutine  we  got  the  3 x 33  polar  tensor  matrix,  F , which 
was  a juxtaposition  of  the  3x3  APTs  of  Che  atoms  of  Che 
molecule  under  investigation! 


9Pjj/3x^  8px/3yi  3pj^/3Sj  ap^^/JXj 


?x  - spy/axj 


*Pv''9/l  *P„/ 


SPj/9Xj  3pj/3yj  9pj/SZj  }p^/ 


where  the  first  three  columns  of  P pertained  to  atom  1,  the 
next  three  to  atom  2,  and  so  forth.  Each  of  the  components 
of  Pjj  were  calculated  in  units  of  e by 

3Pjj/9x.  ■ flp^/dx^  ■)  Pj((Xq  + 0.D05)-  Pjj(Xq  - 0.005))  * 


(4.60324/(2(0.005) )) 


(2-3), 


forth.  Thus 


the  derivatives 


dipole  nonent  with  respect  to  cartesian  di splacenents  of  the 
atoms.  The  cartesian  force  field.  F (converted  to  internal 
symnetcy  coordinates)  and  Che  AFTs,  P , of  2-pyriinidincl  are 
given  in  Chapter  3, 

2. 1.1. 3 RHF  SCF  STO-3G  calculations 

The  subroutine  FREQ  of  GAUSSIAN  82^^  was  used  to 
determine  the  cartesian  force  fields  of  both  tautomers  of 
2-oxopyrimidine . The  calculations  were  done  at  Che  RHF  SCF 
level  with  an  ST0-3G  basis  set.  The  default  parameters  of 
the  subroutine  were  used.  The  reference  geometry  was  the 
optimised  ST0-3G  geometry  of  Section  2. 1.1.1,  This 
subroutine  evaluated  Che  entire  force  field  by  analytic 
gradient  techniques.  No  arts  could  be  determined  with  this 
subroutine.  The  only  necessary  input  was  Che  reference 
state.  It  was  possible  to  use  this  subroutine  because  the 
basis  set  utilised  was  small  (ST0-3G).  This  subroutine  could 
not  be  used  with  bigger  basis  sets  such  as  4-31G  due  to 
computer  limitations.  This  was  the  reason  why  subroutine 
FORCE  was  used  to  evaluate  the  4-31G  force  field  of 
2-pyrimidinol  (Section  2. 1.1. 2}  instead  of  subroutine  FREQ. 
2.1.2  GAUSSIAN  76  Calculations 

Diagonal  force  constants  and  atomic  polar  Censors  for 
various  reference  states  (Chapters  4 - 8)  of  both  tautomers 

were  calculated  at  the  RFH  SCF  level  with  4-31G  basis  sets 
using  the  GAUSSIAN  76  program. The  reference  states  of 


opciinized  geometries 


2-pyrinidinol  were  the  4-31G  and  STO-3G 
(Section  2.1.1.1J.  and  two  possible  crystal  structures 
{Chapter  5).  The  reference  states  of  2-pyrimidinone  were  the 
STO-36  optiaiiaed  geometry  (Section  2, 1,1.1)  and  the  crystal 
structure  (Chapter  6).  The  default  parameters  of  the 
GAUSSIAN  76^^  program  were  used.  The  convergence  was 
5 • 10”^,  and  the  integral  threshold  was  10””  at  all  times. 
The  only  exception  was  the  determination  of  the  diagonal 
force  constants  and  AFTs  of  2-pyrimidinol  at  a convergence 
of  10  ^ (reference  state  is  4-31G  optimized  geometry  — see 
Chapter  5) . 

For  a given  reference  state  the  energies  and  dipole 
moment  components  of  the  ground  state  and  displaced 
geometries  were  evaluated.  From  the  energies  the  force 
constants  were  determined,  and  icon  the  dipole  moments,  the 
AFTs.  It  should  be  mentioned  that  all  calculations  done  with 
GAUSSIAN  76^^  were  carried  out  in  H coordinates.  This  is  a 
set  of  coordinates  described  by  Chln,^^  who  wrote  a secies 
of  programs  in  these  coordinates  that  could  be  used  in 
conjunction  with  GAUSSIAN  76.^^  In  the  present  study  we  call 
them  W coordinates  although  Chin^^  used  a different 
nomenclature.  These  coordinates  are  the  mass  weighted, 
molecule  fixed,  symmetry  adapted,  principal  cartesian 
coordinates.^^  They  are  based  on  a molecule  fixed  coordinate 
system  which  translates  and  rotates  along  with  the  molecule. 
In  this  way  the  vibrational  motions  of  the  molecule  are 
separated  from  translational  and  rotational  ones 


{ Eckact-Sayvetz  condi tions ) . For  no  other  najor  reason 
than  sltnpliclty  in  cooputations  at  that  time  we  used  these 
coordinates . 

For  a given  reference  state  displacements  in  the 
positive  and  negative  directions  were  made  for  each 
coordinate,  and  the  energy  and  dipole  moment  components  were 
caiculated,  one  point  at  a tine,  for  the  reference  state  and 
for  each  displaced  geometry.  The  diagonal  force  constants  of 
the  force  field  in  W coordinates,  F , were  determined  by 


^wi  * • '^'*0  * '^‘*0  - “’■>  -2E(Xg|  I • 

f 4.3598/(0.011^  1 


(2-4) 


where  the  ith  diagonal  force  constant  i of  the  force 

field  was  given  in  ndyn/A  using  displacements  Awi  of  0.01 


displaced  geometries, 
reference  state  in  the 
converted  to  cartesian 
suitable  transformations 

coordinates  in  a man 
displacing  by  Awi)  to 
converted  to  cartesian 


of  the  positively  and  negatively 
and  I was  the  energy  of  the 

sane  units.  The  F matrix  could  be 
coordinates  (X  coordinates)  by 
(see  Section  2.2.1). 
derivatives  were  evaluated  in  H 
analogous  to  Equation  2-3  (but 
the  P matrix.  The  P matrix  was 

■■  '’x'  ’’y 


whec«  P was  in  units  of  e,  u was  the  coocdlnate 
transformation  matrlx^^  (between  H and  X coordinates),  and  H 
was  a diagonal  matrix  containing  the  masses  of  the  atoms  in 
amu.  In  this  fashion  the  diagonal  elements  of  r . and  the 
APTS,  P , could  be  obtained  for  each  reference  state.  The 
results  of  these  calculations  are  shown  in  Chapters  5 and  6. 

2 . 2 Vibrational  Spectra 

This  part  deals  with  the  use  of  ^ Initio  calculations 
of  force  constants  and  APTs  in  the  evaluation  of 
frequencies,  normal  coordinates,  and  intensities.  All  the 

Chin’s  programs. This  study  is  divided  in  four  different 
sections!  interconversion  of  coordinates,  method  to 
determine  frequencies,  determination  of  intensities,  and 
interpretation. 

2.2.1  Interconversion  Of  Coordinates 

As  proposed  by  Pulay  and  coworkers^ ^ we  decided  to 
use  symmetrized  internal  coordinates,  S coordinates,  for  the 
determination  of  IP  frequencies  from  force  constants.  The  S 
coordinates  were  obtained  by  applying  projection  operator 
techniques'^  to  get  linear  combinations  of  internal 
coordinates  corresponding  to  the  Irreducible  representations 
of  the  symmetry  group  of  the  molecule.  All  the  molecules 
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considered  in  the  present  work  have  C syninetry.  The 
corresponding  5 coordinates  were  defined  according  to  the 
recommendations  by  Pulay  and  collaborators.^^ 

The  ab  initio  calculated  force  constants  and  force 
fields  were  expressed  either  in  w or  X (cartesian) 
coordinates.  Since  the  calculations  covered  in  the  following 
sections  were  done  in  S coordinates,  it  was  Important  to 
convert  force  fields  in  w and  X coordinates,  F and  into 
force  fields  in  S coordinates,  F , These  interconversions 
were  carried  out  by  the  following  relationships  (In  matrix 
notation) s 


(2-S) 

(2-7) 

(2-8) 


where  U was  the  coordinate  transformation  matrix  mentioned 
in  Equation  2-5  and  was  its  transpose,  n was  the  diagonal 
matrix  containing  Che  masses  of  Che  atoms  in  amu,  and  A and 


{ between 


B were  coordinate  transfornation  natrices  B 
X coordinates),  and  and  B^  were  their  corresponding 

transposed  natrices. 

The  A and  B matrices  are  quite  important 
relationships.  Their  construction  has  been  explained  by 
several  authors . To  deternine  them  the  first  step 
carried  out  was  to  calculate  the  B.  matrix  which  transforns 
X {cartesian)  to  R (internal)  coordinates: 


R . Bj  * X (2-12) 

where  B^  was  a (3N  a)  * 3N  matrix,  N was  the  number  of 
atoms  of  the  molecule,  and  a was  any  whole  number  (number  of 
redundancies  plus  six).  Thus  6,  could  be  redundant.  To 
eliminate  such  redundancy  and  utilise  the  symmetry  of  Che 
molecule  Che  following  was  used; 


where  H was  a (3N-6)  * (3N  *■  a)  matrix  which  symmetrized  the 
internal  coordinates,  R,  to  give  the  corresponding  S 
coordinates.  In  a similar  fashion  Che  5.  matrix  was 
transformed  by 

B - H • Bj  (2-14) 
where  B transformed  X into  S coordinates  (not  redundant). 


clearly  by  ccnbining  Equations 


through  2-14 


. B • X (2-151. 

The  A natrix  is  a form  of  inverse  matrix  of  B because  it  is 
a rectangular  matrix^^  (A*B«E  and  B*A*S,  where  E is  the 
identity  matrix).  The  following  relationship  was  used 
(further  explanations  on  the  A matrix  can  be  found  in 
references  57  and  59): 


the  inverse  of  the  kinetic  energy  natrix  G.  The  G natrix  is 
defined  in  the  imnediate  next  section. 

2.2.2  method  to  determine  Frequencies 

In  accordance  with  Wilson  et  al.^^  the  G matrix  was 
defined  as 


Furthermore 


(2-171. 


where  F was  Che  force  field  (in  S coordinates  in  our  case), 
L was  the  transformation  matrix  which  related  S coordinates 


diagonal  natcix 


to  normal  coocdinataa  and  A waa  a 
eigenvalues  X.  which  were  related 
frequencies  v.  by 


normal  vibration 


(2-19) 


where  c was  the  speed  of 
particular  wavenumber  in  ccr 
related  to  the  X.  eigenvalue  c 


ht  In  en/sec,  Vj  was  a 
for  the  fundamental  band 
e diagonal  A matrix. 


bands  could  be  obtained  by  the  solution  of  Equation  2-lS.  In 
our  calculations  we  followed  the  solution  of  the  GF  matrix 
of  Equation  2-18  as  explained  by  niyazawa  which  included 


symmetry  considerations . These  methods  have 
various  autors^  and  have  been  detailed  by  See 
2.2.3  Determination  of  intensities 

To  determine  the  intensities  of  bands  t 
formulated  by  Person  and  Newton  were  used^^  ( 


been  worked  by 
t-Lebron . 


i Intensity, 


of 


given  nondegenerate  b 


where  I 3p/9Q^ | was  the  absolute  value  of  the  dipole  moment 
derivative  with  respect  to  the  ith  normal  coordinate. 


ap/SQj  by 


(2-22) 


where  A and  L were  as  previously  defined.  The  was  a 3 ■ 
(3N-6)  matrix  (for  non  linear  molecule  with  N atoms) 
containing  the  dipole  moment  components  with  respect  to 
normal  coordinates  (1  through  3H-€) 


9p,/3Ql  ^Px'^’OsN-e 

8Py/S0i  SPy/SOi  *Py''»03M-6 

3pj/8Qj  8Pj/30,  8pz/9Q3j,_g 


Equation  2-5  could  be  used  to  convert  any  set  of  APTs  in  M 
coordinates,  P , to  P . Once  the  matrix  P was  obtained  it 
could  be  used  to  calculate  P^,  which  in  turn  gave  Che 
intensities  of  the  bands  in  km/mole  via  Equation  2-21. 

2.2.4  Interpretation 

In  the  present  worli  the  normal  coordinates  for 
fundamental  modes  were  explained  in  S coordinates.  These 
coordinates  have  helped  us  to  visualize  the  normal  vibration 
corresponding  to  a given  fundamental  band. 


analyze 


cesults.  The 


first  one  was  derived  by  Horlno 
referred  to  as  the  "MK  method,"  Th 
by  Ksresstury  and  Jalsovsky^^  and  i 
method . " 


and  Kuchltsu®® 
5 referred  to  as 


The  "HR  method"  rationale  was  as  follows,  for  the 
eigenvalue  X|^  of  the  kth  normal  coordinate, 


‘’ik 


Since  the  potential  energy,  V,  was  given  (In  matrix 
notation)  by 


V - O.S  • 0*^  • A • Q (2-25), 

It  followed  that  the  potential  energy  distribution  of  the 
kth  normal  node,  V(0t.)  was 

V(0,)  - 0.5O^^£  Fjj  * ■ L.^  (2-26). 

Horlno  and  Kuchltsu^^  assumed  that  the  product 
^ij^lk^jk  significant  only  when  i-j.  Thus  only  the 

terns  were  evaluated  determine  Che  potential  energy 


dlBtcibutlans  (PEDs).  The  percentage  measuring  the 
contribution  to  the  hth  normal  coordinate  from  the  ith 
ayrametr/  coordinate  was  given  by  PEDjj^, 


where  was  the  diagonal  force  constant  associated  with 
the  ith  symmetry  cordinate,  and  A.  was  the  eigenvalue  of  the 
kth  normal  coordinate.  The  rationale  used  in  the  "KJ  method” 
was  the  same  with  the  difference  that  no  assumptions  were 
made  as  in  the  ”MR  method”;  i.e..  all  terms  in  Equation  2~26 
were  considered.  Thus  the  ”KJ  method"  is  more  theoretically 
justified  than  the  ”HK  method.” 

The  FEDS  give  information  about  the  contribution  of 
each  symmetry  coordinate  to  a given  normal  mode.  In  this 
fashion  one  can  determine  which  5 coordinate  contributes  to 
a given  band.  In  our  calculations  contributions  less  than 
10%  were  not  reported.  Both  methods  have  given  similar 
results  for  the  in-plane  vibrations  of  the  tautomers  of 
2-oxopyrimidine . The  ”RJ  method"  seems  to  bs  mors  accurate 
for  out-of-plane  modes  where  the  assumptions  made  In  the 
”HK  method”  plus  the  low  eigenvalues  could  lead  to 


complications . 


CHAPTER  a 

INTERPRETATION  OF  THE  IR  SPECTRUM  OF  2-PFRIMIDlNOL 

This  chapter  Is  concerned  with  the  Interpretation  o£ 
the  IR  spectrum  of  2-oxopyrimldine  In  a low  temperature 
Argon  matrix.  The  experimental  and  theoretical  techniques 
connected  with  this  spectrum  ace  outlined  first.  Following 
this  the  experimental  fundamental  bands  are  assigned  by 
comparison  with  the  theoretical  calculations  of  the 
spectrum.  Consequently  the  chapter  is  divided  into  four 
sections:  experimental  techniques,  theoretical  techniques, 
assignment  of  bands,  and  summary. 

3 . 1 Experimental  Techniques 

The  experimental  work  was  carried  out  by  M.  J.  Nowak, 
H.  Bostkowska,  and  K.  Sscsepaniak.  A paper  on  these  results 
is  in  preparation.^  Details  of  the  technique  have  been 
previously  outlined . ^ A simple  summary  is  presented  in 


The  2-oxopyrimidine  vapor 
condensed  on  a Csl  window  cooled 
gas  to  solute  ratio  was  about 
(frequencies  and  absorbances)  was 
S80B  spectrophotometer  for  the 
cm  ^ resolution. 


and  the  matrix  gas  were 
to  about  10°  R.  The  matrix 
1000:1.  The  IR  spectrum 
recorded  on  a Per)(in  Elmer 


4000-400 


Table  1-1  lists  Che  experimental  freguenci 


integrated  relative  absorbances  (R.A.s]  obtained  fcon  the 
above  mentioned  spectrum.  The  R.A.s  were  calculated  by 
integration  of  the  bands  in  the  experimental  spectrum  with  a 
compensating  polar  planimeter  of  fixed  arm  (Keuffel  and 
csser  Co.,  Model  * 6200001.  These  R.A.s  are  expressed  in 
unite  such  that  the  integrated  absorbance  of  Che 
experimental  fundamental  vl  was  assumed  to  be  the  same  as 
the  intensity  calculated  for  vl  (•  110  kn/mclel.  The  upper 
parts  of  Figures  1-1  and  3-2  represent  Che  experimental 
spectrum  tabulated  in  Table  1-1.  The  R.A.s  are  drawn  as 
sticbbars  to  facilitate  future  comparisons  with  theoretical 
results.  This  "sticbbar  representation"  is  used  throughout 
the  rest  of  this  dissertation.  As  mentioned  before 
(Introduction),  this  spectrum  is  Icnown  to  be  due  to 
absorption  predominantly  by  the  2-pyrimidinol  tautomer  with 
some  trace  absorption  by  the  2-pyrimidinone  tautomer.^ 

3 . 2 Theoretical  Techniques 

Calculations  were  carried  out  for  2-pyrimidinol.  The 
geometry  was  optimised  with  GAUSSIAN  62^^  using  a 4-31G 
basis  set  (see  Chapter  2).  Using  this  optimized  geometry  as 
reference  state,  the  force  constants  and  apts  were  evaluated 
with  GAUSSIAN  82^^  using  the  4-31G  basis  set  (Chapter  2). 
From  these  results  the  frequencies,  PCDs,  and  intensities  of 
the  fundamental  bands  were  determined  using  the  methods 


outlined 


Section  2.2. 


EXPERIMENTAL  RESULTS  AND  ASSIGNMENT  OF  THE  IR  SPECTRUM 
MEASURED  IN  A DILUTE  AR  MATRIX  FOR  2 -OXOPTRIMIDINE  AT  10°  K 


3S84 

3574 

3053 


2997 


166B 


1627 

1600.5 

1590 


49.2 


V2  12S9,75S10,12S11 

v3  S2S9.47S11 

^J4  35S9.24S10,40S11 

"Region  a" 

"Region  a" 


■Region  b" 


21S2.10S3.23S5, 14 


1573 

1469.5 
1456 
1445 

1434.5 

1390 

1351.5 
1329 


1215 

1211 

1200 

1186 

1168 

1089 


"Region 

"Region 

v8 

"Region 


17S1,11S3,28S4.13S17 
12S2 , loss, 21S7 .15S17 , 12S1B 


1,27S1S,2SS19 


v9  16S3,14S7,40sn 

NlO  22S16.11S17, 36519 

vll  45S16.16S1B 


vll 

vll 

vl2  32S3.34S6.20S1B 

vl3,vl4  19S2.49S5.12S1S  (\13) 


Table  3-2  9ive&  Che  4-31G  optiiriiaed  geometry  of 

2- pycimldlnoI  In  terns  of  coordinates  defined  In  figure 

3- 3.  The  4-316  force  field  in  S coordinates,  P , is  shown  in 
Table  3-3.  The  S coordinates  used  for  3-pyrimidinol  are 
given  in  Table  3-4.  Table  3-5  shows  the  AFTs  in  accordance 
to  the  orientation  indicated  in  Table  3-6. 

Table  3-7  gives  the  unsealed  wavenumbers  icro  and 
intensities  (km/mole)  predicted  for  the  fundamental  bands 
from  the  above  mentioned  force  field  and  set  of  APTs.  The 
PEDs  are  given  in  Table  3-S  for  both  PED  methods  described 
in  Chapter  2.  The  bottom  parts  of  Figures  3-1  and  3-2  show  a 
Stickbar  plot  of  the  unsealed  predicted  spectrum,  to  be 
compared  with  the  experimental  spectrum. 

3 . 3 Assignment  of  Bands 

In  order  to  compare  the  experimental  and  the 
theoretical  spectra  to  reach  an  assignment  of  bands,  it  is 
useful  to  scale  the  calculated  wavenumbers.  This  can  be  done 
because  it  is  known  that  predicted  wavenumbers  ace  too  high 
by  about  104  (see  Introduction).  Consequently  the  calculated 
spectrum  is  modified  slightly  to  obtain  Figures  3-4  and  3-5. 
The  experimental  results  shown  in  these  figures  are 
identical  to  those  shown  in  Figures  3-1  and  3-2.  The 
modification  is  in  the  calculated  spectrum.  All  the 

predicted  wavenumbers  are  multiplied  by  a constant  factor  of 
0.9  (Table  3-7),  leaving  the  Intensities  untouched. 

This  modification  of  predicted  wavenumbers  is  done 
only  for  assignment  purposes.  It  has  proven  to  he  quite 
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ZED  GEOHETBY  OF  2-PYBIHIDINOL 
AT  TRE  HF  SCF  LEVEL  HITH  A 4-31G  BASIS  SET 


10 

11 


0.9S3 

1.068 

1.069 


14 

15 

16 


21 

22 

24 

25 

26 


117.2 

122.1 

116.4 

121.4 

112.4 


121.9 

116.5 


116.2 


122.1 


*8ond  lengths,  c,  and  Angles.  B,  ace  aceocdlng  to  Figure  3-3. 
The  geometry  was  optimised  by  the  BFGS  procedure  of  Head  and 


Stcucture,  3tom  numbering  and  in-plane 
internal  coordinates  of  2-pyrimidinol  used 
in  this  work.  The  out-of-plane  internal 
coordinates  are  the  out-of-plane  bends  and 
torsions  given  in  Table  3-i, 


so 


51 


TABLE  3-3 

FORCE  CONSTANTS  FOR  2-PTRIMIDINOL 
IN  SYMMETRIZED  INTERNAL  COORDINATES  S 
BVALOATED  AT  THE  HF  SCF  LEVEL  WITH  A 4-31G  BASIS  SET 
g^a,b,e,d  gj  gj  g^  gj  gg 


6.90261 

1.29878 

-.75255 

.72709 


6.92856 

.95023 

-.75839 


95745 

84444 

06873 

01695 

02918 

02536 

54593 

17644 

65215 

07127 

03076 

03879 

03426 


75171 

90136 

05266 

01912 

03137 

01439 

46613 

37366 

54017 

06073 

04073 

03936 

02910 


8.83356 

1.14173 

-.66902 

1.011Q9 


-.01774 

-.02835 

.23403 

.01695 

-.56756 

-.00376 

-.36224 

.00276 

-.00492 


7.93078 
.90268 
-.68362 
-.06565 
.01163 
. 08635 
.06409 
-.00846 
.00660 
.06268 
.37299 
.08953 
.02506 
.16025 

-.01656 


-.15368 
-.00189 
-.00801 
. 06532 
.08560 


8.72396 

-.04195 

-.00117 

-.02663 

.19956 


14079 

07933 

01831 

01762 

18557 

15622 


46933 

32839 

08317 

00751 

00316 

00027 

35361 


SI  7.91SZ8 
sa  .03642  8.93739 


0007S 

01605 

23673 

24604 

37339 

11380 

49751 


S17 


-.00971 


SIB  -.00343 
S19  .00839 


-.00098 

-.00276 

-.00808 

-.01362 

.02281 

.08132 

.29046 

-.00064 

-.00100 

-.00259 


01368  6.38005 
00229  .01357 
08529  -.11521 
08673  .05342 
02147  -.09266 
00918  .00140 
00167  -.00215 
01156  -.00448 
00162  .00067 
00771  .00398 


6.31713 

.08425  1.70178 
.02429  .01515 

.08548  -.00648 
-.00888  -.02262 
-.00182  .02196 

.00749  -.01526 
-.00184  -.00117 
-.01263  .01456 


S19 


513  1.67077 

514  -.16276 

515  -.06862 

516  .00651 

517  .04143 

518  -.09422 

519  .08831 


1.83607 

-.00657  1.37391 

-.00108  .16156 

.07999  -.01220 

-.05470  -.00692 

-.00521  -.01425 


.89260 

.00322  .69191 
.00224  .01008 
.00219  -.01193 


.60824 

.00988  .68520 


cominued 


.34371 


824  .00630  -.00177  -.01384  -.02474 

S25-. 15007  .16175  .02646  -.00342 


S4 

TABLE  3-4 

SYHHETBIZED  INTERNAL  COORDINATES  FOR  2-FYRIHIDINOL 


SS  • AcS 


8S 


tC4C5 

rC5C6 

rC6Nl 

rC20 


SIO 


arlO 


a/2)f. 

(1/2)“ 

U/2)“ 

(1/2)® 

(1/2)® 


’(a9i2-asi3+ 

•®(2flS12-4pi 


'®(as20-fie24i 

■®{i021-i025) 

■®|a022-fl026) 


s0i4-a0i54^A0ie-a0iT) 

3-a0i4+2fl0i5-fl0i6-aani 

6-4017) 


iC20 

JCOH 

6C4H 

iCSH 


TftBLB  3-4  continued 


S2C  - U/6l®'^l6Tl-4T2*AT3-flT4*AT5-flt6) 

521  - (l/2)|atl-4T3<-a-t4-flT6) 

522  - (l/12)®'^(-4Tl-f2aT2-4T3-at4t24T5-4T6) 

523  - aroa 

524  - (1/2)® '^f 4x7+418) 

525  - 4vH8 
S2€  - 4vH10 
S27  - avHll 


iH9 

rmi 


^The  nomenclature  and  numbering  follow  Figure  3-3. 

The  recommendations  of  Pulay  et  al.  (Reference  11)  are 

^The  atoms  involved  in  torsions  are:  ilw6123,  i2wl234r 
T3-234S,  T4-34S6,  iS-4561,  and  TS-SS12. 


TABLE  3-S 

ATOMIC  POLAR  TENSORS  POR  2-FYRlMlOINOL 
EVALUATED  AT  THE  HF  SCF  LEVEL  HITH  A 4-31C  BASIS  SET 


-.4S981  -.12B3S  0 

.36156  -1.03260  0 

0 0 -.40389 


.23131  -.07765  0 

-.02949  .49910  0 

0 0 .45602 


-.00369  -.05791  0 


°The  aton  numbering  is  given  in  Figure  3-3;  the  coordinate  system 
orientation  is  the  principal  cartesian  coordinates  of  Table  3-6. 

‘^The  reference  geometry  is  the  4-31S  optimized  geometry  of 
2-pyrimidinol  (Table  3-2). 
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helpful  in  einiler  s i tuatione . Thus  the  dssignnent 
can  be  carried  out  using  Figures  3~A  and  3-S. 

The  fundanental  bands  in  Che  experlDental  spectruin  are 
split  in  various  instances,  e.g.,  fundamental  v5  (Figures 
3-4,  3-5).  We  shall  not  attempt  to  define  the  exact  causes 
for  this  experimental  splitting.  In  general,  band  splitting 
is  due  to  matrix  effects,  hydrogen  bonding,  Fermi  resonance, 
tunneling,  and  rotational  isomerism.  Further  complications 
are  found  when  tautoraerism  and  polymeric  aggregation  are 
taking  place  significantly.  As  mentioned  in  Section  1.1  many 
of  these  effects  can  be  discarded  as  far  as  this  dilute  Ar 
matrix  spectrum  is  concerned.  The  2-oxopyrimidine  system  is 
known  to  be  only  in  the  monomeric  2-pyrimidinol  form.^'^  One 
can  safely  rule  out  tautomerism,  hydrogen  bonding  and 
polymeric  aggregation,^'^  Hence  Fermi  resonance,  rotational 
isomerism,  and  matrix  effects  are  considered  Co  be  major 
reasons  for  Che  experimental  band  splitting  observed  here. 

In  the  experimental  spectrum  of  2-oxopyrimidine  small 
bands  that  could  not  be  easily  explained  have  been  grouped 
together  as  "region  a,"  "region  b,"  "region  c,"  and  "region 
d"  respectively  (see  Figures  3-1,  3-2,  3-3,  and  3-4). 

Possible  explanations  for  the  absorption  in  these  regions 
are  combinations  and/or  overtones  of  fundamental  bands  of 
2-pyrimidinol,  the  presence  of  traces  of  2-pyrimidinone , and 
possibly  impurities. 

we  now  concentrate  on  specific  regions  of  Che 
experimental  spectrum  and  compare  Che  absorption  there 
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theoretical  predictions 


Co  make  the  assignments. 
The  PEDs  by  the  "KJ  method"  (Table  3-fl)  are  used  to  describe 
the  normal  nodes  of  the  fundamental  bands. 

3.3.1  aeqlon  <000-2000  cm~^ 

Two  bands  around  3600  cm  ^ and  a series  of  weak  bands 


around  3100-3000  CBi~  are  seen  experimentally  in  this  region 
(see  Figure  3-4).  Based  on  our  prediction  for  vl  the  two 
experimentally  observed  bands  around  3600  cm  ^ are  assigned 
to  a split  vl  0-R  stretching  node.  Possible  reasons  for  the 
splitting  are  matrix  effects,  tunneling,  and  rotational 
isomerism.  The  first  three  experimental  weak  bands  observed 
around  3100-3000  cra~^  are  assigned  to  the  C-H  stretching 
modes  vZ,  v3,  and  v4  respectively,  node  v2  is  due  to  a 
mixture  of  all  three  C-H  stretches  where  rCSK  is 

predominant.  Mode  v3  is  a Isl  combination  of  C4H  and  C6H 
stretches.  Mode  v4  corresponds  to  a roughly  equal 

combination  of  all  in-plane  C-H  stretches.  The  bands  in 
"region  a"  (3013  - cm“^)  are  probably  due  to 

combination  bands  involving  v6  and  v6  modes  which  have 
strong  absorptions.  These  results  are  consistent  with  the 
assignments  based  on  characteristic  frequencies.^^  in  this 
region  the  assignments  for  other  purine  and  pyrimidine 
bases^  ^ are  similar  to  Che  ones  proposed  here.  So  other 
fundamentals  are  expected  in  this  region  for  2-pyrimidinol . 
3.3.2  aeqion  2000-14D0  cm  ^ 

In  this  region  (Figure  3-S)  the  experimental 

absorption  bands  assigned  to  the  vS  through  v6  fundamentals 
are  observed.  These  bands  are  commonly  referred  to  as 


breaching  vibracione  of  aromatic  cinge.^^  The  weak 
experimental  bands  of  "region  b"  (1722.5  - 160D.S  cm 

"region  e"  [H58  - 1445  cm"^).  and  "region  d"  (1417  - 1390 

cm  are  also  seen  here.  By  comparing  the  upper  and  lower 
parts  of  figure  3-5  it  can  be  seen  that  the  assignment  of 
experimental  bands  is  straightforward.  The  only  complication 
concerns  Che  u5  mode.  The  calculated  intensities  for  the  uS 
and  u6  fundamentals  are  roughly  equal.  Thus  the  experimental 
doublet  assigned  to  v5  is  believed  Co  be  due  to  splitting  of 
the  band  by  Fermi  resonance  so  that  the  experimental 
intensity  for  u5  is  the  sum  of  the  intensities  of  bands  in 
the  doublet.  The  ratio  between  the  calculated  intensities 
for  v7  and  v8  is  opposite  to  Che  one  observed  for  the 
corresponding  experimental  bands.  This  fact  can  be 

attributed  to  the  limitations  of  the  level  of  theory  and 
basis  set  used  in  the  calculation.  It  is  possible  that 
because  of  these  limitations  the  aromatic  electronic 
delocalization  is  not  properly  explained.  This  in  turn  would 
affect  the  values  of  force  constants  and  intensities.  The 
FED  results  (Table  3-8)  indicate  Chat  modes  v5  through  v8 
are  due  to  In-plane  bends  and  ring  stretches.  Bands  in 
"region  b"  (Figure  3-5)  are  probably  due  to  C*0  stretches  of 
2-pyrimidinone  tautomer  and/or  Co  combination  bands  of 
2-pyrimidinol . ^ Similarly  bands  in  "region  c"  and  "region 
d"  are  probably  combination  bands^^  of  lower  frequency  modes 
of  2-pyrimldinol . 


vacioue  other  pyridine  and  pyrimidine  aystema  show 
similar  breathing  modes  and  combination  bands  in  this 
region. A point  of  interest  is  the  split  v5  band.  Since 
many  combination  bands  are  possible  near  this  band,  rermi 
resonance  seems  a plausible  reason  for  the  splitting.  The 
latter  explanation  is  preferred  over  matrix  effect  or 
rotational  isomerism  arguments.  This  is  inferred  from  the 
fact  that  the  normal  coordinates  (PEDs)  for  all  breathing 
modes  are  similar  and  all  should  be  split  for  the  same 

A series  of  experimental  absorption  bands  are  observed 
in  this  region  (Figure  3-5).  From  comparison  of  calculated 
spectra  with  observed,  the  assignment  of  v5  and  vlO 
fundamentals  ie  straightforward.  These  bands  ace  due  mostly 
to  couplings  of  in-plane  C-H  bends  with  the  C-0  stretch  and 
to  a lesser  extent  with  the  COH  bend.  The  vll  and  vl2 
fundamentals  are  harder  to  assign,  He  believe  that  vll 
fundamental  in  the  experimental  spectrum  is  split  into  five 
bands.  The  reason  is  because  the  intensity  calculated  for 
vll  is  too  high  to  be  explained  by  one  of  the  low  intensity 
bands  observed  in  this  pact  of  the  spectrum.  The  FED  table 
shows  that  vll  is  predominantly  the  in-plane  COH  bending. 
This  band  nay  be  expected  to  be  quite  sensitive  to  the 
environment  of  the  molecule.  The  vll  band  has  no 
contribution  from  in-plane  fCOH,  and  is  expected  to  be  a 
single  band  quite  in  agreement  with  our  calculations. 


Four  weak  exparinenCal  bands  are  observed  around  1100 
cm  ^ (Figure  3-SI.  Calculations  in  this  area  predict  two 
in-plane  fundamentals,  vl3  and  vl4,  and  two  very  weak 
out-of-plane  fundamentals,  v20  and  v21.  The  corresponding 
experimental  bands  are  so  weak  and  close  to  each  other  that 
we  assign  then  to  the  vl3  and  vl4  fundamentals.  Based  on  the 
oaloulated  intensities  for  the  out-of-plane  fundamentals  v20 
and  v21  we  do  not  expect  Chat  they  be  seen  at  all  in  Che 
experimental  spectrum.  The  FED  table  shows  that  vl3  and  vl4 
are  due  mostly  to  in-plane  ring  stretches,  and  v20  and  v21 
to  out-of-plane  H bends,  tH’s. 

The  splitting  of  the  experimental  vll  fundamental  is 
of  particular  interest.  To  explain  it  in  terms  of  Fermi 
resonance  or  combination  bands  does  not  seem  right.  This  is 
evident  by  looking  at  the  experimental  out-of-plane  v2f 
fundamental.  This  band  is  due  mostly  to  out-of-plane  torsion 
involving  the  0 and  H atoms  of  the  hydroxyl  bond,  tOH,  The 
absorption  assigned  to  v23  displays  a very  similar  pattern 
to  that  for  vll  band  which  is  due  to  in-plane  4COH.  Since 
Fermi  resonance  or  combination  bands  (not  possible  at  all  in 
the  v2S  area)  cannot  explain  the  structure  for  the  v25  band 
it  could  be  argued  that  they  also  do  not  explain  either  the 
analogous  vll  band.  This  leaves  us  with  matrix  interaction, 
tunneling,  and  rotational  isomerism  (i.e.,  in-plane  or 
out-of-plane  cotamers  of  the  OH  group)  as  the  most  probable 
reasons  for  the  splitting  of  the  vll  band.  The  OH  group  can 
be  visualised  as  subject  to  various  orientations  with 
respect  to  the  ring,  and/or  interactions  within  the  matrix. 


3.3.4  Region  1000-200  cn~^ 

In  this  region  e few  experinentel  low  frequency 
in-plane  fundamentals  and  most  of  the  experimental 
out-of-plane  fundamentals  are  observed  (Figure  3-5).  A weak 
band  predicted  for  vl5  around  1000  cm~^  due  Co  in-plane  ring 
bends,  S6I,  is  not  seen  experimentally.  In  Che  900-8D0  cm~^ 
portion  of  the  spectrum  the  experimental  bands  are  assigned 
as  ul5.  u22,  and  u23  respectively  based  on  Che  corresponding 
predicted  bands  in  this  region.  Halogenation  studies,^  which 
replaced  HIO  of  2-pyrinidinol  by  Br  or  Cl,  showed  that  the 
experimental  fl7I  cm’^  band  of  the  spectrum  of  2-pyrinidinol 
remains  unchanged  throughout  halogenation.  This  was  not  the 
case  for  the  810  cra~^  band.  Prom  the  FED  results  (Table  3-8) 
\22  is  shown  Co  depend  more  on  the  coordinates  of  the 
molecule  in  the  region  immediately  adjacent  to  HIO;  hence  it 
is  sensible  to  assign  it  to  the  band  observed  experimentally 


The  assignment  of  Che  remaining  experimental 
fundamentals  is  rather  straightforward  by  comparison  with 
the  calculated  spectrum.  The  major  difficulty  occurs  for  the 
v25  fundamental  due  to  the  out-of-plane  torsion  tOH.  This 
split  band  has  been  already  explained  in  connection  with 
discussion  of  the  vll  fundamental  involving  th  in-plane 
bend  of  the  COH  group.  The  predicted  v27  band  cannot  be 
compared  to  experimental  results,  since  the  experimental 
study  did  not  cover  the  region  below  400  cm 


Table  3-1  sumraari aes  the  assignment  and  Intecpretation 
of  the  IR  spectcun  of  2-oxopycinidine  in  a dilute  low 
temperature  Argon  matrix.  For  purposes  of  simplifying  these 
results  Table  3-9  is  constructed.  This  table  provides  a 
direct  comparison  of  experimental  to  calculated  spectra  for 
2-pyrimidlnol , The  experimental  frequencies  are  the  assigned 
fundamentals,  disregarding  peaks  from  regions  a,  b.  c.  and 
d.  For  split  bands  such  as  vS  the  average  of  the  split  band 
frequencies  is  listed  as  the  expe cinental  value.  For  vl3  and 
vl4  the  Same  frequency  is  reported,  the  experimental 
absorbances  are  the  R.A.s  reported  in  Table  3-1.  For  split 
bands  the  R.A.s  are  added  so  as  to  give  the  total  R.A.  for 
the  fundamental.  For  vl3  and  vl4  the  R.A.  values  are  assumed 

The  calculated  frequencies  reported  in  Table  3-9  are 
the  calculated  scaled  frequencies  {multiplied  by  the  single 
scaling  constant  of  0,9)  as  described  in  Table  3-7.  The 
calculated  R.A.s  of  Table  3-9  are  the  unsealed  calculated 
intensities  (Table  3-7)  multiplied  by  a conversion  constant. 
He  call  this  conversion  constant  the  absorbance-intensity 
conversion  constant.  The  absorbance-intensity  conversion 
constant  for  the  in-plane  case  is  the  ratio  of  the  sum  of 
in-plane  experimental  R.A.s  to  Che  sum  of  the  in-plane 
calculated  intensities  (Table  3-7).  It  is  equal  Co  0.S77 
(•772.1/1336.0).  Similarly  the  absorbance-intensity 
conversion  constant  for  Che  out-of-plane  modes  is  0.360 
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(•121.2/336.6).  The  unobserved  v27  is  not  taken  Into 
account.  This  conversion  constant  is  essentially  the  product 
of  the  cell  path  length,  the  solute  concentration  and  a 
scaling  constant  averaged  over  the  number  of  considered 
fundamentals.  Thus  the  usage  of  this  conversion  constant  can 
hardly  be  called  scaling  of  the  absorbances.  It  is  used  to 
enable  us  to  have  a more  proper  comparison  of  the  band 
absorbances.  This  method  can  only  be  used  after  the 
assignment  of  the  experimental  spectrum  has  been  made. 
Figure  3-6  shows  the  plot  of  the  results  of  Table  3-9. 

The  agreement  shown  in  Figure  3-6  and  Table  3-9  show 
that  the  relationship  between  experimental  and  calculated 
spectra  is  quite  satisfactory.  The  only  significant 
modification  of  the  calculated  results  is  the  multiplication 
of  all  frequencies  by  a single  scaling  constant,  0.9.  The 
discrepancies  that  were  seen  between  the  R.A.s  and  the 
intensities  of  figures  3-4  and  3-9  are  minimised  by  the 
usage  of  the  absorbance-intensity  conversion  constants.  The 
methods  used  in  the  construction  of  Table  3-9  will  be  useful 
in  future  chapters. 

It  is  interesting  that  the  absorbance-intensity 
conversion  constant  used  for  the  out-of-plane  modes  is 
smaller  than  for  the  in-plane  modes.  This  is  consistent  with 
previous  findings  in  calculations  of  intensities  of  large 
molecules  at  the  HP  6CF  level  with  small  basis  sets.  It  was 
found  that  in  general  out-of-plane  intensities  are  predicted 
to  be  much  higher  in  comparison  to  experimental  values  than 
are  those  predicted  for  the  in-plane  intensities. 
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In  short  it  can  be  said  that  the  assignnent  ie 
satisfactory,  and  in  good  agreenent  with  literature  results 
throughout  the  entire  IR  region.  The  usage  of  the  4-31G 


basis  set  at  the  bhf  scf  level  is  rewarding,  It  shows  that 
this  level  of  theory  can  be  used  reliably  for  large 
molecules  like  2-pyrioidinol . Figures  3-1  and  3-2  which 
compare  the  experimental  frequencies  and  intensities  to  the 
unsealed  calculated  spectrum  of  2-py r imidinoi  could  have 
been  just  as  well  used  for  assignment  purposes  as  were 
Pigures  3-4  and  3-5.  These  figures  are  included  in  this 


level  quite  satisfactory  unsealed  results  are  obtained.  The 
usage  of  0.9  as  a scaling  constant  for  frequencies  ie  done 
exclusively  Co  ease  the  work  of  assignment  of  experimental 
spectra.  This  is  quite  different  from  procedures  to  fit 
calculated  Co  experimental  spectra  by  statistical  methods. 


figure  3-6.  Coioparison  ot  experimental  to  calculated 
spectra  of  2-pyrinidinol . The  5XPT  values 
correspond  to  analysis  of  experimental 
results  (Table  3-91.  The  CALC  values 
correspond  to  analysis  of  full  4-31C 
study  (Table  3-9 } . 
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EFFECT  OF  UPGRADING  THE  FORCE  FIELD  IN  DIFFERENT  COORDINATE 
SYSTEMS  ON  PREDICTED  IR  SPECTRUM  OF  2-PYRIHIDINOl 

Ab  discussed  in  Chapter  1 it  is  customary  to  upgrade 
or  scale  force  fields.  No  consistent  study  has  been  carried 
out  on  the  effect  on  predicted  IR  spectra  of  choosing 
different  coordinate  systems  for  the  scaling  of  force 
constants.  The  purpose  of  this  chapter  is  to  study  this 
effect.  The  presentation  will  be  divided  in  three  ma^or 
sections:  methods,  survey  of  results,  and  conclusions. 

4 .1  Methods 

The  optimised  geometry  of  2-pyrimidinol  calculated  at 
the  RF  SCF  level  with  an  ST0-3G  basis  set  (Section  2. 1.1.1) 
was  utilised  as  reference  state  to  evaluate  the  force  field 
for  this  molecule  at  the  HF  SCF  level  with  an  STO-3G  basis 
set  (Section  2.1.1.31.  Table  4-1  gives  this  optimized 
geometry.  Table  4-2  shows  this  force  field  expressed  in 
symmetrized  internal  coordinates  (S  coordinates). 

This  8TO-3G  force  field  was  upgraded  by  scaling 
according  to  Equation  1-3  where  was  this  ST0-3G  force 

field;  C was  the  diagonal  matrix  of  scaling  constants  chosen 
to  be  the  ratios  of  the  diagonal  force  constants  of  the 
4-31G  force  field  (Table  3-3)  to  the  corresponding  values  of 
the  STO-3G  force  field  (Table  4-2).  For  example  for  SI  the 
conversion  constant  was  a . 90281/0 . B657S  - 1.D042.  These 

conversion  constants,  also  referred  to  as  scaling  constants. 


OPTIHIZED  GEOHETHY  OF  2-PYRIHIDINOL 
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METHODS 


THE  4-31G/STO-3G  SCALING 
FOX  2-PYXIHlDINOL  USED  IN  THE  SC, 
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and  the  force  field,  F was  obtained.  Force  field  F Is 
referred  to  as  a "near  4-llG"  force  field. 

This  "near  4-31G"  force  field  was  then  used  to 
calculate  band  frequencies  and  potential  energy  distribution 
{PEDs)  of  1-pyrlnidinol  by  the  method  to  determine 
frequencies  and  the  "HK  method"  respectively  of  Chapter  2. 
The  band  intensities  were  calculated  via  Equation  2-22  where 
the  L matrix  from  this  "near  4-31G"  force  field  was  used 
together  with  the  atomic  polar  tensors  (APTs)  of 
2-pyrimidinol  evaluated  at  the  HF  5CF  level  with  a 4-31G 
basis  set  (Table  3-S).  This  whole  process  is  defined  by  the 
name  "SC  method"  which  stands  for  the  "symmetrised  internal 
coordinate  method."  The  frequencies,  PEDs,  and  intensities 
obtained  by  this  method  are  given  in  Tables  4-4,  4-5,  and 
4-6  respectively  under  the  SC  heading. 

Similarly  we  considered  other  prescriptions  for 
upgrading  the  5TO-3G  force  field  expressed  in  other 
coordinate  systems.  The  "HC  method"  and  the  "XC  method"  were 
done  with  an  F^^  force  field  (Equation  1-3)  calculated  for  w 
and  X coordinate  displacements  respectively.  The  W 
coordinates  were  the  ones  described  in  Chapter  2,  and  the  X 
coordinates  were  the  cartesian  coordinates.  The  F^^  force 
field  used  with  the  HC  and  XC  methods  was  the  same  STO-3G 
force  field  of  Table  4-2  expressed  in  w and  x displacement 
coordinates  respectively.  These  coordinate  interconversions 
are  shown  in  Equations  2-6  through  2-11. 
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■thods  of  scaling  ace  desccibed  in  Section 


The  4-31G  di&9onel  values  used  to  evaluate  the 
4-31G/STO-3G  scaling  constants  vets  obtained  fcoiri  the 
diagonal  values  of  the  4-31G  force  field  of  2-pyrimidinol 
(Table  3-3)  expressed  in  W or  X coordinates.  These  4-31G 
diagonal  values  In  H or  X coordinates  were  divided  by  Che 
corresponding  diagonal  values  of  the  ST0-3G  force  field 
(Table  4-2),  F**'  of  Equation  1-3,  in  w or  X coordinates 
respectively.  Thus  the  4-31G/STO-3G  scaling  constants  which 
made  up  the  diagonal  C matrix  of  Equation  1-3  were  obtained. 
In  this  way  Che  "near  4-31G”  force  field,  F of  Egustlon  1-3, 
in  W or  X coordinates  was  obtained.  In  order  to  simplify  the 
discussion,  the  conversion  constants  for  the  wc  and  xc 
methods  are  not  Included  in  the  present  woric.  They  can  be 
calculated  from  available  data  by  Equations  2-6  through 
2-11.  using  the  "near  4-31G"  force  field  in  (4  or  X 
coordinates  and  the  4-31G  APTs  (Table  3-S ) the  band 
frequencies,  PEDs,  and  intensities  of  2-pyrimidinol  were 
calculated  in  the  same  way  as  for  the  "SC  method,”  and  these 
values  are  included  in  Tables  4-4  through  4-6. 

It  should  be  realized  Chat  although  scaling  in  the  SC, 
WC,  and  XC  methods  was  done  via  Equation  1-3,  the  force 
field  to  be  upgraded  was  a matrix  with  a different  degree  of 
diagonality  in  each  case.  In  order  Co  eliminate  (to  some 
extent)  this  differing  factor,  two  other  methods  were 
deviced,  the  "mCd  method"  and  the  "XCD  method." 

In  the  case  of  the  "MCD  method"  the  "near  4-31G"  force 
field  obtained  from  the  "WC  method"  was  converted  from  W 
into  S coordinates  (via 


Equations  2-6  through  2-11).  Then 


the  diagonal  values  of  this  "converted  near  411-G"  force 
field  were  divided  by  the  corresponding  diagonal  values  of 
the  STO-3G  force  field  in  s coordinates  {Table  4-2J  to  give 
4-31G/STO-3G  scaling  constants.  These  constants  are  referred 
to  as  and  ace  given  in  Table  4-3.  using  Equation  1-3 
with  the  values  defining  the  C natrix.  the  STO-3G  force 
field  in  S coordinates  (Table  4-2)  for  the  "near  4-31G" 
force  field,  F,  in  S coordinates  corresponding  to  the  "WCD 
Dethod"  was  obtained,  using  this  "near  4-31G"  force  field 
and  the  4-31G  APTs  (Table  3-S)  the  band  fcequencles,  PEDs, 
and  intensities  of  2-pyrimidinol  were  calculated  in  the  same 

For  the  "XCP  method"  the  prescriptions  were  identical  except 
that  X coordinates  were  used.  The  values  of  Table  4-3 
correspond  to  Che  4-31G/STO-3G  scaling  constants  used  in 
this  method.  The  scaling  constants  and  results  from  the  WCD, 
XCD,  and  SC  methods  can  be  directly  compared  since  they 
involve  the  upgrading  of  the  identical  f'"^  force  field 
expressed  in  5 coordinates,  i.e.,  they  upgrade  a force  field 
with  the  same  degree  of  diagonality. 

Tables  4-4  and  4-5  show  the  band  frequencies  and  PEDs 
of  2-pyrimldinol  obtained  from  Che  usage  of  the  "near  4-31G" 
force  fields  for  the  various  methods.  These  tables  include 
two  extra  headings,  STO-3G  and  4-31G,  which  are  added  for 
comparison  purposes.  The  ST0-3C  results  were  obtained  by 
using  Che  unsealed  STO-36  force  field  of  2-pyrimidinol 
(Table  4-2)  to  determine  band  frequencies  and  PEDs  by  the 


method  to  determine  frequencies  and  the  "HK  method'  oC 
Section  2.2.  The  4-31G  results  are  taken  directly  from 
Tables  3-7  and  3-8,  Table  shows  the  band  intensities 

calculated  for  2-pyrimidinol  by  the  various  methods  using 
the  4-31G  APTs  from  2-pyrimldlnol  {Table  3-3). 

4 . 2 Survey  of  Results 

Figures  4-1  through  4-S  show  the  spectra  predicted  by 
the  various  methods  considered  in  this  chapter.  These 
spectra  are  stickbar  plots  of  Tables  4-4  and  4-6.  The 
breakdown  Into  the  various  spectral  regions  is  arbitrary. 
These  figures  demonstrate  that  the  various  upgrading  methods 
are  not  equally  effective  In  imitating  the  4-31G  spectrum. 
For  instance  in  Figure  4-1  it  Is  deduced  by  visual 
Inspection  that,  although  all  methods  improve  on  the  STO-3G 
results,  the  "SC  method'*  seems  to  give  Che  best  agreement 
with  4-31G  spectrum.  Comparison  of  the  performance  of  Che 
various  methods  by  visual  inspection  of  the  spectra  tends 
not  to  be  very  objective.  In  order  to  achieve  more 

objectivity  in  the  comparison  of  the  ability  of  the  methods 
to  imitate  4-31G  spectrum  simple  statistical  approaches  are 
considered . 

Table  4-7  contains  a simple  statistical  analysis  of 
the  efficiency  of  the  various  methods  In  reproducing  the 
4-31G  results.  The  efficiency  is  measured  by  the  success  in 
mimicking  the  band  frequencies,  PBDs,  and  intensities  of 
2-pyrimidinol  evaluated  at  the  HF  SCF  level  with  a 4-31G 
basis  set,  i.e.,  4-31G  results,  The  EF,  EN,  and  Ei  terms 


Cooiparieon  of  calculated  IR  spectra 
through  v4)  from  various  scaling 
methods.  The  methods  are  elaborated 
Chapter  4. 


92 


93 


STATISTICAL  ANALYSIS  OH  THE  ABILITY  OF  THE  VARIOUS  METHODS 
TO  MIMIC  THE  FBEOICTEQ  4-31G  RESULTS  FOR  THE  vl  THROUGH  v4 
FUNDAMENTAL  BANDS  OF  2-FYRIHlDINOL 
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full  4-31G  study  of  Chapter  3. 


frequency,  norpal  ncdee,  and 


neaeuce  deviations  in 
intensities  respectively  from  the  4-31G  results.  For 

instance  EFi  is  the  difference  in  frequency  between  the  vl 
band  predicted  by  a qiven  method  and  the  vl  band  predicted 
from  the  4'31G  results  (Table  4-4).  The  lEF,  £EN,  and  CEI 
terms  are  Che  sun  of  the  absolute  values  of  the  deviations 
of  predicted  frequencies  (EFs),  normal  modes  (CNs),  and 
intensities  (Els)  respectively  for  a given  method  for 
fundamentals  vl  through  v4.  The  FP  is  an  assigned  frequency 
parameter  which  helps  one  to  judge  Che  methods  from  the  best 
{•  1)  to  Che  worst  (-6)  in  their  ability  to  mimic  the  4-31G 
results  on  the  basis  of  EEF  values.  Similarly  N?  and  IP  are 
the  normal  mode  and  intensity  parameters  based  on  grading 
the  IBM  and  £BI  values  respectively.  The  CFF  is  the 
efficiency  parameter  which  is  used  to  grade  Che  overall 
performance  of  a given  method  based  on  the  FP,  NP,  and  IP 

Table  4-7  permits  us  to  make  a more  objective  and 
quantitative  comparison  of  spectra  of  Figure  4-1  obtained  by 
different  methods  Chan  just  a visual  comparison.  Similarly 
Tables  4-6  through  4-11  permit  us  to  compare  the  performance 
of  each  of  the  methods  for  the  spectra  calculated  in  Figures 
4-2  through  4-5  respectively.  Table  4-12  summarises  these 
tables  where  the  analysis  is  divided  according  to  the 
symmetry  of  the  molecule  (in-plane  and  out-of-plane 
fundamental  modes).  In  this  table  the  EFIP,  EMIP,  and  BIIP 
are  the  means  of  the  lEF,  £EN,  and  ESI  values  respectively 
for  the  in-plane  fundamentals  (vl  through  vl9).  The  EFOP, 


Comparison  of  calculated  is  spectra  (vS 
through  uB)  from  various  scaling 
methods.  The  methods  ace  elaborated  in 
Chapter  4. 
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Conpaclson  of  calculated  IR  spectra  (v9 
through  ul5)  Scorn  various  scaling 
fflethods.  The  nethods  are  elaborated  in 
Chapter  4 . 


Conparison  of  calculated  IR  spectra 
through  vl9)  fron  various  scaling 
methods.  The  methods  are  elaborated  In 
Chapter  A . 


102 


103 


Figure  4-S . 


Comparison  of  calculated  IR  spectca  (u20 
through  »27)  from  various  scaling 
methods.  The  methods  ere  elaborated  in 
Chapter  4. 
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values  respectively  tor  the  out-of-plane  fundamentals  (v20 
through  v27).  The  EPTOT,  ENTOT,  and  SITOT  are  the  means  of 
the  lEF,  £EN,  and  EEl  values  respectively  for  all 
fundamentals  (vl  through  v27).  The  FPIP,  NPIP,  IPIP,  and 
SFFIP  are  the  frequency,  normal  node,  intensity,  and  overall 
parameters  analogous  to  Tables  4-7  through  4-11  which  help 
one  to  grade  the  performance  of  the  various  models  regarding 
In-plane  modes  only.  Similarly  PPOP,  NPOP,  IPOP,  and  EFFOP 

out-of-plane  modes  only;  and  FPTOT,  NPTOT,  IPTOT,  and  EPFTOT 
are  used  to  grade  the  performance  of  the  models  throughout 
all  fundamental  nodes. 


4 . i conclusions 

The  various  methods  considered  in  this  chapter  gave 
predicted  results  of  different  quality  for  band  frequencies, 
normal  modes,  and  intensities  of  2-pyrinidinol . The  reason 
for  these  differences  did  stem  from  the  use  of  different 
coordinate  systems  in  upgrading  procedures.  To  actually 
judge  as  to  how  severe  these  differences  are  is  hard  to 
determine  from  the  data  presently  available  to  us.  The 
interesting  tact  is  that  the  same  force  fields  at  the  hf  SCF 
level  with  4-nC  and  ST0-3G  basis  sets  respectively  were 
used  in  all  the  methods  considered  in  this  chapter.  This 
finding  alone  shows  us  that  there  is  the  need  to  carry  out 
further  investigations  on  the  selection  of  coordinate 
systems  to  be  used  in  scaling  procedures. 


Hithin  the  scope  o£  the  statistical  approaches  used  so 
far  the  trend  in  the  deviations  of  the  nocnal  nodes  and 
intensities  seens  to  be  sinilar  to  the  trend  in  the 
deviations  in  the  band  frequencies  throughout  Che  various 
nethods  applied  Co  calcuiations  for  the  in-plane  inodes. 
These  results  are  not  entirely  without  contradictions 
especially  when  one  exanines  the  out-of-plane  results. 
Nonetheless  it  seems  certain  that  the  conunonly  used  approach 
to  consider  only  the  comparison  of  band  frequencies  in  the 
evaluation  of  Che  results  of  different  scaling  procedures 
for  force  fields  is  insufficient.  It  is  possible  that  more 
sophisticated  statistical  analysis  is  needed  to  investigate 
Che  differences  between  the  predicted  spectra  of  the  various 
methods.  Ideally  this  analysis  should  look  simultaneously  at 
the  differences  in  band  frequencies,  normal  modes,  and 
intensi ties . 

The  analysis  in  Table  4-12  for  the  in-plane  results 
indicates  that  the  best  method  is  SC.  For  the  out-of-plane 
results  XC  is  Che  preferred  method.  All  upgrading  methods 
show  an  improvement  over  the  "ST0-3G  method”  which  uses  an 
unsealed  force  field  at  the  HF  SCF  level  with  an  STO-3G 
basis  set.  By  judging  the  efficiency  over  all  the 
fundamental  modes  (EFFTOTI  it  seems  that  SC  is  the  method  of 

Table  4-12  contains  a few  other  parameters  that  have 
not  been  explained  yet.  The  HDIP  and  MDOP  parameters  are  the 
means  of  the  absolute  magnitudes  of  the  differences  between 


the  and  values  from  the  values  of  Table  4-3  for 
the  in-plane  and  out-of-plane  modes  respectively.  The  HIP 
and  HOP  pacametecs  are  the  means  of  these  differences  for 
in-plane  and  out-of-plane  nodes  respectively.  The  3D1AG1P 
and  %DIAGOP  are  cough  measures  of  the  degree  of  diagonality 
of  the  in-plane  and  out-of-plane  portion  of  the  force  field 
to  be  upgraded,  i.e..  the  force  field  of  Equation  1-3. 

The  (DIAGIP  is  the  ratio  of  the  sum  of  the  absolute  values 
of  the  diagonal  elenente  of  the  in-plane  force  field  to 
the  sum  of  the  absolute  values  of  all  the  in-plane  elements 
of  the  force  field,  multiplied  tines  100.  The  tblAGOP  is  the 
same  corresponding  ratio  for  the  out-of-plane  force 

constants. 

The  purpose  of  introducing  these  parameters  is  to  show 
that  there  is  some  relationship  between  then  and  all  the 
rest  of  statistical  parameters  concerning  frequencies  and 
normal  nodes  in  Table  4-12.  The  fact  Chat  XC  and  HC  methods 
had  values  of  the  tDIAGIP  AND  tDIAGOP  quite  different  from 
Chose  for  the  SC  and  STO-3G  methods  was  an  important  factor 
in  Che  developnent  of  the  WCD  and  XCD  nethods.  Thus  the  XC 
and  WC  methods  are  not  included  in  these  immediate 

observations.  On  the  other  hand  Che  SC.  XCD,  and  HCD  nethods 
share  the  same  values  of  Che  3D1AG1P  and  tDIAGOP.  There 
seens  to  be  a relationship  between  Che  errors  for  the 
predicted  in-plane  frequencies  and  for  Che  normal  nodes 
predicted  by  these  nethods  and  the  corresponding  HDIP  and 
HIP  values.  This  sane  trend  is  harder  to  follow  in  the 
out-of-plane  case.  These  observations  could  be  roughly 
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interpreted  ee  eaylng  that  the  closer  the  scaling  constants 
of  a given  method  ace  to  the  C values  the  better  the  method 
will  mimic  the  4-31C  band  frequencies  and  normal  modes. 
Notice  that  the  band  intensities  are  not  included  In  this 
discussion  because  they  depend  not  only  on  the  force  field 
but  also  on  the  APTs  (Section  1.2.5). 

Various  observations  are  made  in  this  chapter.  The 
need  exists  to  understand  scaling  methods  more  thoroughly. 
It  seems  that  the  more  diagonal  the  force  field  to  be 
upgraded  is,  the  better  the  results  Chat  are  obtained  from 
scaling  procedures.  Thus  s search  could  be  suggested  for 
coordinate  systems  that  provide  a high  degree  of  diagonality 
for  the  force  field  to  be  scaled.  The  usage  of  other  systems 
of  conversion  could  be  tried  with  perhaps  the  inclusion  of 
the  geometry  of  the  molecule.  As  seen  before  band 
Intensities  and  normal  nodes  should  be  Incorporated  In 
statistical  fitting  of  calculated  Co  experimental  spectra. 

The  philosophy  used  in  this  chapter  la  to  try  to 
imitate  the  spectrum  of  2-pyrimidinol  calculated  at  the  HP 
SCF  level  with  a 4-51G  basis  set  using  the  results  at  the 
sane  level  calculated  with  an  5T0-3G  basis  set.  It  is  shown 
that  such  imitation  is  not  easily  achieved  by  the  various 
methods  attempted.  Further  still  the  relationships  between 
these  two  levels  of  theory  is  not  clearly  understood  as  to 
their  effect  on  predicted  Ih  spectra.  If  this  is  true  for 
this  comparison  between  two  rather  similar  results 
calculated  with  rather  similar  representations  of  the 
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waveCunctionB . the  situation  should  be  even  mote  complex 
when  a theoretical  representation  is  scaled  to  fit 
experimental  results.  It  is  suggested  that  such  conparisons 
between  calculations  from  rather  similar  levels  of  theory 
te.g.,  using  4-31G  and  ST0-3G  basis  sets  as  here},  be 
studied  further  In  order  to  understand  the  theory  of  scaling 
of  force  fields. 


jrroundlng  various  types 
work  as  discussed  in 
some  of  our  findings  in 
le  predicted  IB  results 


CHAPTER  S 

EFFECT  OF  CONVERGENCE,  INTEGRAL  THRESHOLD,  AND  CHOICE 
OF  REFERENCE  STATE  IN  AB  INITIO  CALCOLATIONS  ON 
THE  PREDICTED  IR  SPECTBOM  OF  2-PYRIHIDINOL 

There  are  nany  controversies  s 
of  calculations  in  predictive  IP 
Chapter  1.  This  chapter  presents 
connection  with  these  problems.  1 
from  ab  Initio  calculations  perforraed  at  different  degrees 
of  convergence  and  integral  threshold  are  discussed. 
Similarly  the  effect  of  choosing  different  reference 

polar  tensors  are  considered.  All  the  studies  covered  in 
this  chapter  are  made  exclusively  on  2-pyrinidinol . The 
discussion  is  divided  in  three  sections:  convergence  and 
integral  threshold,  reference  geometry,  and  conclusions. 

5.1  Convergence  and  Integral  Threshold 
Ab  initio  calculations  of  diagonal  force  constants  and 
atomic  polar  tensors  of  2-pyrinidinol  were  carried  out  at 
the  HF  SCF  level  with  a 4-31G  basis  set  at  two  different 


! limits  of  the  GAUSSIAN 
• 10  The  integral  threshold  i 
Diagonal  force  constants  and  f 
convergence  are  referred  to  as  HC 
atomic  polar  tensors  respectively, 
ID  ^ convergence,  they  are  referred 
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,s53  , 


10'®  and  5 


program: 
i both  cases  was  10'® 
’Ts  evaluated  at  lo' 
force  constants  and  H 
Similarly  for  the  5 


and  atonic  polar  tensors  respectively.  The  details  of  the  ab 
initio  calculations  are  given  in  Section  2.1.2.  The 
reference  geonetry  was  at  all  times  the  4-31G  optimized 
geometry  of  2-pyriinidinol  (Tables  3-2  and  3-6). 

By  the  '*HCD  method"  explained  in  Chapter  4.  the  HC  and 
LC  diagonal  force  constants  and  the  ST0-3G  force  field  of 

2- pyrimidlnol  (Table  4-2)  were  used  to  give  HC  and  LC  "near 
4-31G"  force  fields  respectively.  The  scaling  constants  are 
given  in  Table  S-1,  where  C^^  and  C^^  correspond  to  the  HC 
and  LC  cases  respectively.  The  HC  and  LC  atonic  polar 
tensors  of  2-pyrimidinol  are  given  in  Tables  S-2  and  5-3 
respectively. 

Using  the  HC  and  LC  "near  4-31G"  force  fields,  the 
band  frequencies  of  Tables  5-4.  5-5.  and  the  normal  nodes  of 
Table  5-B  were  obtained.  The  4-31G  APTS  (Table  3-5),  the  HC 
APTs,  and  the  LC  APTs  of  2-pyrinidinol  were  used  together 
with  HC  and  LC  "near  4-31G"  force  fields  respectively  to 
give  the  band  intensities  shown  in  Tables  5-4  and  5-5. 
Tables  S-6  and  S-7  are  added  for  comparison  purposes.  The 
band  frequencies  in  Table  5-6  are  the  unsealed  frequencies 
obtained  from  the  4-31G  force  field  of  2-pyrimldlnol  (Table 

3- 3).  The  band  intensities  of  Table  5-6  are  obtained  fron 
this  4-31G  force  field  and  the  4-31G  APIs  (Table  3-5),  the 
HC  APTs.  and  the  LC  APTs  respectively.  Similarly  Table  5-7 
shows  the  band  frequencies  and  intensities  obtained  from  the 
unsealed  STO-3G  force  field  of  2-pyr imldinol  (Table  4-2), 
and  the  above  mentioned  APTs.  In  short,  four  different  force 
fields  (HC,  LC,  4-31G,  and  STO-3G)  were  used  in  the 
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evaluation  of  frequencies  and  pegs.  For  each  force  field 

three  sets  of  APTs  H-31Gr  HC,  and  LC)  were  used  in  the 

evaluation  of  intensities. 

Table  S-9  sunnarises  the  statistical  analysis 
perforned  on  these  results.  The  errors  in  frequencies,  PEOs, 
and  intensities  are  calculated  in  the  same  way  as  presented 
In  Chapter  4.  These  errors  ace  measured  from  the  deviations 
of  the  various  cases  from  4-31G  results  for 
2-pyrimidinol  (Chapter  3).  Since  three  sets  of  APTs  are 
being  considered,  three  measures  of  EIIP  and  ElOP  ace  given. 
The  HDIP,  MOP,  MDIP,  and  KOOP  values  are  measured  with 
respect  to  c scaling  constants  (Table  4-3}  as  was  done  in 

Chapter  4.  These  values  can  be  compared  to  each  other  since 

the  force  field  to  be  upgraded  had  the  same  diagonality  in 
all  cases.  The  statistical  results  for  the  "SC  method" 
(Table  4-12)  ace  added  for  comparison  purposes. 

Table  S-9  is  a useful  summary  that  shows  the 

differences  in  prediction  between  the  various  considered 
methods  and  the  4-31G  results  of  Chapter  3.  Loo)ring  at  the 
frequency  results  foe  HC  and  LC  one  can  deduce  that  the 
higher  the  convergence  the  better  the  results.  Nevertheless 
the  actual  difference  in  the  EFIP  for  these  two  cases  is  not 
as  large  as  the  corresponding  difference  in  EFOP.  This 
suggests  that  convergence  is  more  important  for  the 

determination  of  the  out-of-plane  force  field  than  for  the 
in-plane  one.  At  least  it  can  be  said  that  the  the  effect  of 
convergence  is  important  for  the  out-of-plane  modes  where 
the  use  of  C^^  scaling  constants  leads  to  a set  of  predicted 
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TABLE  S-4 

PREDICTED  BAND  FREQUENCIES  AND  INTENSITIES 
FOR  IR  FUNDAHEHTAL  BANDS  OF  2-PYRII1IDINOL  FROM 
HC  'NEAR  4-31G"  FORCE  FIELD  AND 
4-aiG,  KC,  AND  LC  ATOMIC  POLAR  TENSORS 


3950.4 

34Si!2 

3439.2 

1032!4 

isiela 

1346]e 

1084i? 

1011. a 

73li7 

457!o 


1336.8 
667. S 


26 

27 


1181.3 

1166.9 

963.6 

927.2 


101.8 

10.3 

6.1 


107.9 

12.4 


0 

119.1 

9.7 

203]6 

16.8 


356.3 


lEl'  59.3 


51.1  54.4 


®Fcequencle8  wete  predicted  fcon  HC  "near  4-31G’  force 
by  the  "HCO  method"  (Section  5.1). 

field^and^a  glven'set'of'*APTs^°4-31o*APTs  (Table  1-5), 
HC  APTS  (Table  5-2),  and  LC  APTs  (Table  5-3). 


field 


TABLE  5-5 

PREDICTED  BAND  FRSCUENCIE6  AND  INTENSITIES 
FOR  IB  FUNDAMENTAL  BANDS  OF  2-PYBIHIDINOL  FROM 
LC  "NEAR  4-310"  FORCE  FIELD  AND 
4-31G,  KC,  AND  LC  ATOMIC  POLAR  TENSORS 


1 3965. 7 

2 3491.0 

3 3462.9 

5 190713 

6 1875.2 

7 1730.3 

9 IS34I9 

10  ISIS. 3 

11  1367.0 

12  1262.7 

13  1234.9 

14  1112.3 

15  1040.0 

16  973.0 

17  7S2.6 

16  637. B 

19  472.4 


142.1 

6.3 

18.5 

I66I7 
247.9 
179.6 
192. S 


25I7 

23l3 

1337.8 


20  1244.4 

21  1219.7 

22  1027.8 

23  989.1 

24  715.9 

25  684.8 

26  544.4 

27  286.3 


112.6 

.4 

121.8 

s!4 

4.8 

336.8 


142.8 

6.1 

18.1 

166!? 

248.0 

179.4 

192.3 


0 

121.2 

.2 

117.5 

88.9 

338i6 


1366!i 


‘Frequencies  are  predicted  Cron  LC  "near  4-31G"  force  Held 

“intensities  ace  predicted  by  usage  of  LC  "neat  4-31G 
field  and  a given  set  of  APTs:  4-31G  APts  (Table  3-5 


r-1, 


TABLE  5-6 

PREDICTED  POTENTIAL  ENERGT  DISTRIBUTION 
FOB  ZR  FUNDAMENTAL  BANDS  OF  2-PIRiniDZNOL  FROM 
"NEAR  4-31G"  FORCE  FIELDS  AT  TWO  CONVERGENCES 


10S9 

10S9 

22S1 

1251 


2SS4.1 
2657 

12S7,24Sia,14S19 

1587,21517,17819 

sssia 

15517,12819 

69812 


10S17 

1181. 

1084, 

S2S16 

1487, 

16816 


1385,2787 
1186, 1187, 23S18,1 
22S17,16S19 
21817,17819 
■■  ■■  27819 


2282,2687,12814 

72513,17814 

16S7,10S13,S9S14 


20  64825,43826,21827 

21  37325,78827 

22  36820,35823,70826 

23  126520,45823,29825,31827 

24  17820,15321, 39S22,6BS23, 
25826 

25  93824 

26  101521,27522,10825 

27  24821,79822 


100S2 

67820 

36827 


29826 

14826,95527 

1,65823,46526 

15823,31825,29826, 


40823,37324,17326 


104821,30822,11825 
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frequencies  Chan  Is  worse  Chan  that  obtained  utilisinq  an 
unsealed  ST0-3G  force  field.  This  is  a form  of  warning  that 
caution  should  be  exercised  in  scaling  procedures,  where  in 
some  cases  scaling  a force  field  actually  gives  worst 
results  Chan  not  scaling  at  all.  The  difference  in  the  SPIP 
and  EFOP  values  between  Che  SC  and  HC  methods  is  not  large 
(they  differ  mainly  in  integral  threshold  values).  This 
suggests  that  better  integral  thresholds  improved  the 
results  but  not  by  much.  The  trends  concerning  normal  modes 
are  harder  Co  see  since  the  statistical  values  ace  quite 
close  to  each  other.  The  statistical  values  on  the 

KPIP,  HIP,  HOOP,  and  HOP  values. 

The  results  for  predicted  intensities  in  Table  5-9  are 
more  dlfficults  to  follow,  One  can  begin  by  looking  at  the 
APTs  of  Tables  5-2  and  5-3.  These  tables  show  the  EXPO 
parameter  which  stands  for  the  sum  of  the  absolute  values  of 
the  differences  between  the  in-plane  components  of  the  APTs 
and  the  corresponding  components  of  the  4-31G  APTs  of 
2-pyrimidinol  (Table  3-5).  Similarly  COFO  is  the  sum  of 
these  deviations  for  the  out-of-piane  components.  Tt  is 
interesting  Co  notice  there  is  a large  change  in  EIPD  value 
between  the  KC  and  LC  APTs,  whereas  a small  one  is  seen  for 
the  lOPD  case.  The  small  llPD  for  HC  APTs  suggests  that 
integral  threshold  was  not  important  for  the  in-plane  case, 
since  4-31G  and  HC  APTs  had  roughly  the  same  convergence 
(10"®  and  10'®  respectively)  but  quite  different  integral 
thresholds  (10'^®  and  I0~®  respectively).  By  similar 
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reasoning  the  large  CIPD  for  the  LC  A?Ts  (obtained  at  S * 
10  convergence)  suggests  that  convergence  was  quite 
important  for  the  in-plane  case.  For  the  HC  APTs  the 
magnitude  of  the  COPI>  is  similar  to  the  IIPD  perhaps 
suggesting  again  that  integral  threshold  was  not  a 
significant  parameter  in  these  ^ initio  calculations.  The 
fOPD  values  for  HC  and  LC  APTs  are  quite  similar  indicating 
that  convergence  was  not  as  important  as  in  the  in-plane 
case  (at  least  in  Che  10  Co  5 * I0~^  convergence  range). 

Table  5-4  shows  that  the  above  mentioned  obaervatlons 
on  APTs  are  imitated  rather  closely  by  the  predicted  band 
intensities.  The  parameters  £I  and  £EI  in  this  table  are  the 
sum  of  the  intensities  and  the  sun  of  the  absolute  values  of 
the  deviations  of  the  intensities  from  the  corresponding 
4-31G  Intensities  (Table  3-7)  respectively.  For  the  in-plane 
nodes  both  4-31G  and  HC  intensities  were  similar  whereas  LC 
showed  a larger  deviation.  For  Che  out-of-plane  modes  it 
seems  there  was  no  difference,  or  an  element  of  randomness, 
between  the  £BI  values.  The  same  situation  is  observed  in 
Tables  S-5  through  5-7, 

Loo)cing  bacic  at  the  errors  in  intensity  values  of 
Table  5-d,  it  can  be  seen  that  sinilar  observations  to  the 
above  nentioned  ones  are  observed,  For  a given  force  field, 
HC  for  instance,  the  EIIP  and  EIOP  values  follow  Che  same 
trends  throughout  the  range  of  considered  APTs.  The 
logistics  of  understanding  the  statistical  parameters  of  the 
intensities  for  various  force  fields  with  a given  APT  is 
more  difficult  to  develop.  For  the  in-plane  modes  results 
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TABLE  S-9 

STATISTICAL  ANALYSIS  ON  THE  ABILITY  OF  VARIOUS 
FORCE  FIELDS  AND  ATONIC  POLAR  TENSORS  FROK  DIFFERENT 
LEVELS  OF  CONVERGENCE  AND  INTEGRAL  THRESHOLD  TO  HIKIC 
THE  PREDICTED  4-31G  RESULTS  FOR  THE  IR  FUNDAMENTAL  BANDS 
OF  2-PYRIHIDINOL 
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seemed  to  get  better  with  force  fields  of  higher 
convergence.  We  recall  that  band  intensities  were  obtained 
from  the  interplay  of  force  fields  and  APTs  (Eguation  2-22). 
Thus  instances  could  have  occurred  where  the  intensities 
calculated  with  a lower  convergence  force  field  were  similar 
to  the  ones  from  a higher  convergence  force  field.  For 
example  consider  the  rather  close  BIIP  value  for  HC  force 
field  with  HC  APTs,  and  the  value  for  LC  force  field  with  LC 
APTs.  The  reasons  behind  these  observations  on  intensities 
predicted  by  various  methods  need  be  investigated  further. 

All  the  above  mentioned  observations  on  frequencies 
normal  nodes,  and  intensities  predicted  by  the  various 
methods  are  further  summarised  in  the  FP,  NP,  and  IP 
parameters  for  in-plane  and  out-of-plane  nodes,  in  general 
it  can  be  said  that  the  higher  the  convergence  and  the 
integral  threshold  the  better  the  results.  Nonetheless  In 
the  range  of  studied  convergences  and  integral  thresholds, 
the  convergence  effect  seems  to  be  a much  more  critical  one 
than  the  integral  threshold. 

5.2  Befetenee  Geometry 

Ah  initio  calculations  of  diagonal  f 
atomic  polar  tensors  of  2-pyrimidinol  s 


1-31G 


•t  using  three 
program  at  5 * 


reference  geometries  with  the  GAUSSIAN 
10  ^ convergence  and  10  ^ integral  threshold  {see  details 
In  Section  2.1.2).  The  reference  geometries  are  labelled  as 
3GG,  XLl,  and  XL2.  The  force  constants  evaluated  using  these 
reference  geometries  are  referred  to  as  3GG,  XLl,  and  XL2 
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force  constants  respectively,  similarly  APTs  evaluated  with 
these  geometries  are  3GG,  XLl.  and  XL2  atomic  polar  tensors 
respectively.  The  3GG  geometry  is  the  ST0-3G  optimized 
geometry  of  2-pyrimidinol  (Table  4-1).  The  XLl  and  XL2 
geometries  are  possible  crystal  structures  of  2-pyrinidinol . 
The  experimental  geometry  of  2-pyrimidinol  crystal  was  not 
known  because  in  the  crystal  the  2-pyrimidinone  tautomer 
predominates.  The  XI.1  georaetry  was  obtained  by  using 

crystallographic  data  for  pyrimidine^®  together  with  the 
optimized  geometry  from  a calculation  with  an  STO-3G  basis 
set  for  2-pyrimldlnol  given  previously  in  Table  4-1.  The 
ring  portion  was  determined  from  Che  crystallographic  data, 
and  the  exocyclic  components  from  the  theoretical  geometry. 
The  Xt2  geometry  was  determined  by  pursuing  a more  obstinate 
approach  to  estimate  Che  "real"  experimental  geometry  of 
2-pyrimidinol.  The  ring  portion  and  CH  angles  were  obtained 
from  the  pyrimidine  crystal,^®  the  CH  bond  lengths  from 
considerations  on  experimental  work  on  benzene . 
the  CO  and  COH  angles,  the  CO  bond  length  from  crystal  work 
on  6-bromo-2-hydroxypyridi ne , and  the  OH  bond  length  from 
the  considerations  of  Rogers  for  methanol  and  ethanol.®^  The 
XLl  and  XL2  geometries  are  shown  in  Table  5-10. 

By  the  "WCD  method"  explained  in  Chapter  4,  the  3GG, 
XLl,  and  XL2  diagonal  force  constants  and  the  ST0-3G  force 
field  of  2-pyrimidinol  (Table  4-2)  were  used  to  give  3GG. 
XLl,  and  XL2  "neat  4-31G"  force  fields  respectively.  The 
scaling  constants  are  given  in  Table  5-11  where  Cjgg, 
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TWO  POSSIBLE  CRYSTAL  GEOMETRIES  OF  2-PYBIHlDINOL 

XLl  XL2  XLl  XL2 


1 1.335 

2 1.340 

3 1.341 


1.340 

1.341 


126. B 
116.3 


116.3 


16  116.6 


116.6 


18  116.6 

19  103.7 

20  121.3 


11 


1.009 


1.000 


22  116.3 

23  116.3 

24  116. B 


26 
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and  correspond  to  the  3GG,  XLl,  and  XL2  respectively. 

The  AFTe  evaluated  with  these  reference  geometries  are  given 
in  Tables  S-12  through  S-14.  Tables  5-lS  through  5-17  show 
the  corresponding  orientation  used  in  the  calculation  of 

Using  the  3GG,  XLl,  and  XL2  "near  4-31G"  force  fields, 
the  band  frequencies  and  normal  nodes  of  Tables  5-18  and 
S-19  were  obtained.  The  3GG  APTs  were  used  together  with  the 
3GG  "near  4-31G"  force  field  to  give  Che  3GG  intensities  of 
Table  5-18.  Similarly  Che  XLl  APTs  and  XLl  "near  4-31G” 
force  field  gave  the  XLl  intensities  of  Table  5-16;  and  the 
XLl  APTS  and  XL2  "near  4-31G"  force  field  gave  the  XL2 
intensities  of  Table  5-18. 

Table  5-20  summarizes  the  statistical  analysis  of 
these  different  calculations.  The  errors  in  frequencies, 
PEDs,  and  intensities  are  calculated  in  the  same  way  as 
presented  in  Chapter  4.  These  errors  are  the  deviations  of 
the  given  parameter  from  Che  result  for  Chat  parameter 
calculated  with  the  4-31G  basts  set  for  2-pyrlmidinol 
(Chapter  3).  The  MDIP,  HIP,  HDOP,  and  HOP  values  ace 
measured  with  respect  to  C scaling  constants  (Table  4-3)  as 
was  done  in  Chapter  4.  These  values  can  be  compared  Co  each 
other  since  the  force  field  to  be  upgraded  had  the  same 
diagonality  in  all  cases.  The  results  of  the  LC  and  ST0-3G 
methods  with  LC  APTs  from  Table  5-9  are  included  in  this 
table  for  comparison  purposes.  They  had  the  same  convergence 
and  integral  threshold  as  Che  3GG,  XLl,  and  XL2  calculations 
presented  in  Table  5-20. 
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TABLE  S-11 

THE  4-31G/STO-3G  SCALING  CONSTANTS  FOB  2-PTBtMIDlNOL 
USING  VARIOUS  REFERENCE  STATES 


20 

21 

22 

23 

25 

26 
27 


7B78 

0912 

9124 

81S1 

0926 


8346 

0907 

8619 


624S 

6039 

6938 

7470 


IWl 


7779 

6840 

6806 

6974 


7677 

7401 

7621 


1.5620 

1.4203 

1.4206 

2.7268 

l]2151 

1.2197 


1.0397 

.9682 

1.1930 


1.0162 

1.0242 


. 5438 
1.2601 
1.1039 


5421 

6246 

6862 


infocmatlon  see  Section  S^'l. 
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ATOMIC  POLAR  TENSORS 
IBIDINOL  WITH  CRYSTAL  XI 


GEOMETRY 


.18461  .15553  0 

-.04S14  -1.55481  0 

0 0 -.63235 


-.10690  .09361 


.07362  .00456  0 

0 0 .18922 


1.06440  -.36616 

.20552  .67917 


Of  Table  3-6. 
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TABLE  S-14 


"f"  4-.J.33 

^"1-  ,J.. 

-i“  -^1"  ...L 3 

-P” 

■”P"  :#=  ,.3:„ 

-I-  .J33 

^”T’  ^4-  .3,”  . 

-P”  -■!“  -.1, 

-:"■  -"I”  ,J.. 

:li:‘ 

X*  Z® 


• -337.03778011  ( 5 • 10~®  CG,  10'®  IT.  4-31G  BS)'* 


- 2.6563  (10‘®  CG,  10  IT,  STO-3G  BS)' 

. 4.2169  (S  • 10'®  CG,  10'®  IT,  4-31G  BS)'* 


-.274446 

-.835771 

-.199754 

1.140399 

1.833041 

1.064570 

2.215037 

2.479175 


137 


X®  Y®  2® 


- 3.7967  (5  • 10“^  CO,  i0‘®,  4-310  BS ) 
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Table  S-20  helps  one  to  inspect  the  variations  in 
results  frost  the  various  cases  with  respect  to  the  4-31G 
results  of  Chapter  3.  Looking  at  the  frequency  parasieters  in 
this  table  It  can  be  seen  that  the  different  reference 
states  gave  different  deviations.  The  LC  results  were  the 
best  In  general.  This  Is  to  be  expected  since  the  reference 
state  in  this  case  was  the  4-31G  optlnised  geonetry  (Table 
3-2).  The  normal  nodes  for  all  cases  were  approximately 
equivalent.  The  HOIP,  hip,  hdop.  and  HOP  values  mimic  to  a 
good  extent  the  observations  on  frequencies  and  normal 

a set  of  new  parameters  is  Introduced  in  Table  5-20. 
These  are  HDR  and  HR  which  stand  for  the  mean  of  the 
absolute  values  and  the  mean  of  the  values  respectively  of 
the  deviations  of  the  bond  lengths  of  a given  geometry  from 
the  corresponding  bond  lengths  of  the  4-31G  optimized 
geometry  (Table  3-2).  Similarly  HDA  and  HA  correspond  to 
deviations  In  the  angles.  These  are  rough  estimates  of  how 
close  a given  geometry  was  to  the  4-31G  optimized  geometry 
(Table  3-2).  One  night  have  expected  that  the  closer  the 
geometry  of  the  reference  state  was  to  the  4-31G  optimized 
geometry,  the  closer  the  predicted  IR  results  would  have 
been  to  4-31G  results  of  Chapter  3,  i.e.,  better  statistical 
parameters  of  Table  5-20.  It  seems  this  was  not  the  case. 
For  example,  XL2  was  closer  to  the  4-31G  optimized  geometry 
than  3GG  was.  Nonetheless  XL2  deviations  in  force  constants 


TABLE  5-18 

PREDICTED  BAND  PRSOOBNCIES  AND  INTENSITIES  OF  2-PTFIMIDlNOL 
EVALUATED  WITH  DIFFERENT  REFERENCE  STATES 


3699 

3280 

3271 

3255 

1533 


1396.1 

1310.0 

109s!i 

1070.4 

905.9 


20  1334 

21  1270 

23  1052 

24  866 

25  713 

26  575 

27  295 


7 197.5 

4 1.6 

9 3.0 

7 2.0 


993.9 
889.1 

687i6 

451.9 
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TABLE  5-19 

PBEDICTED  POTENTIAL  ENERGY  DISTRIBUTION  FOB  2-PYRIKIDINOL  IR 
FUNDAMENTAL  BANDS  OBTAINED  FROM  DIFFERENT  REFERENCE  STATES 


1 100S6 

2 98S10 

3 96S11 


17S16 

21S18 


IBSl 


14S19 


18S19 

13S6 


25S17 

1SS19 


11  12S1S 
24S17 
26S19 

12  16S6 
53S18 


13 


2SS4 

21S5 

19S17 


16S1 

15S19 


43S10 

20S11 

49S9 

47811 


13S1S 

21S17 

12S19 

16S4 

2355 

12815 


22817 

1SS19 

12S16 

21817 

27819 

1283 

1BS6 

S3S18 

2584 

2385 

1BS17 

10819 


1586 

1587 
15817 
19819 
13817 
1681B 
20819 
12S3 
1385 
2286 
3BS16 
1083 
3284 
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TABLE  5-19  c 
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17S1 

10S3 

67Sl' 

6732 

4353 

60S4 

46SS 

S0S6 


ias2 

25S7 

12S12 


74S13 

19S14 

1SS7 


58S14 
19  92315 


20  109S25 
18526 

21  28526 
96S27 


33S26 

32S20 

23S25 

45S26 

38527 

100524 


25 


22521 

51S22 


26  101S21 


27 


11S25 

70S22 


6751 

6332 

4253 


4655 

4856 


1351 

75512 

7931 

5452 

4653 


4755 

4156 


2057 

10514 


60514 

90515 

97525 

28526 

19526 

93527 


89515 

38526 

77527 

66525 


23527 

138520 

93523 

14525 

16527 

99524 


91523 

16525 

20522 

52523 

36525 

66526 

22520 

20521 

S6S22 

13527 

63520 


40522  20525 

50523  16526 

27526 

105521  99524 

30S22 

12525 

22521  89520 

81522  33521 

24522 
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and  frequencies  were  greater  Chan  3GG  results.  This  suggests 
that  approxlpate  geometries  used  as  reference  states  need 
not  give  approximate  results.  More  is  said  about  this  later. 

The  errors  in  the  intensities  seen  in  Table  S-20  are 
more  difficult  to  understand.  This  is  a consequence  of 
having  used  not  only  different  "near  4'31G”  force  fieids, 
but  also  different  APTs.  As  mentioned  before,  it  is 
difficult  to  evaluate  Che  interplay  between  the  force  field 
and  APTS  that  takes  place  in  the  evaiuacion  of  band 
intensities.  The  best  results  pertained  to  LC  case  as 
expected.  But  this  is  not  necessarily  a fair  comparison.  A 
better  comparison  couid  be  carried  out  by  Che  transfer  of 
APTS  between  reference  states  for  intensity  calculation 
purposes.  This  is  not  to  be  considered  in  the  present  work. 
The  next  section  addresses  this  probiem  in  a different  way. 

5 . 3 Conclusions 

Convergence,  integral  threshold.  and  choice  of 
reference  state  play  a significant  role  in  determining 
predictions  of  IR  spectra.  Cases  are  seen  where  poor 
convergence  or  choice  of  a reference  state  leads  to  results 
that  are  worse  than  unsealed  ones.  i.e..  calculations  with 
an  ST0-3G  basis  set.  It  is  observed  that  in  general  the 
better  Che  convergence  and  integral  threshold  Che  better  the 
results.  The  similarity  between  trends  in  the  EF  and  EN 
values  and  those  in  the  HDIP,  HIP.  HOOP,  and  HOP  values 
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LC 

3GG 

XLl 

XL2 

STO-3G 

:Si% 

v,-.i 

V.i\ 

ll:i 

s;:; 

W.l 

iV.l 

24.4  41.1  S7.9  57.7  37.3 

Isi  ill  i‘i-  il-  1H‘ 

.2808  .5866  .2167  -.3124  .0354 

i i if  if 

; 5 < I ‘ 


Calculations  using  sinilac  rafersnce  states  need  not 


necessarily  predict  similar  IR  spectra.  This  indirectly 
contradicts  researchers  who  use  reference  states  close  to 
the  experimental  geometry  in  the  evaluation  of  force  fields 
expecting  predicted  results  close  to  experimental 
observations  (see  introduction).  On  the  other  hand  this 
indirectly  supports  researchers  who  use  theoretically 
optimized  reference  states  in  the  determination  of  force 
fields  at  the  same  level  of  theory.  Our  conclusions  tend  to 
malce  us  side  with  the  latter. 

k small  observation  should  be  made  at  this  point 
concerning  XL!  calculations.  The  ab  initio  calculation  of 
diagonal  force  constants  in  w coordinates  (at  the  KF  SCF 
level  with  a 4-31G  basis  set)  with  XL2  as  reference  state 

0.05  md/A)  for  two  out-of-plane  nodes.  They  were  set  equal 
to  zero.  This  was  the  only  case  throughout  all  calculations 
carried  out  in  the  present  work  where  negative  force 
constants  were  seen.  This  can  he  interpreted  in  the  sense 
that  at  least  for  fundamentals  with  very  small  frequencies 
it  might  be  dangerous  to  use  reference  states  such  as 
experimental  ones.  It  is  paradoxical  that  the  only  case 
where  this  happened  was  for  XL2,  the  geometry  we  expected  to 
be  closest  to  the  experimental  one. 

Table  5-21  is  constructed  so  as  to  permit  a different 
comparison  between  the  various  cases  treated  in  this 
chapter.  So  far  all  cases  are  compared  to  the  4-31G  results 
of  Chapter  3 during  the  statistical  analysis.  In  this  table 
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the  comparisons  are  made  with  respect  to  expecinental 
results.  The  calculated  frequencies  of  a given  case,  e.g.. 
HC  predictions,  were  multiplied  by  0.9  as  was  done  in 
Chapter  3.  These  scaled  frequencies  were  then  compared  tc 
the  esperlnental  frequencies  of  Table  3-9,  The  CFIF  and  EFOF 
values  were  the  means  of  the  absolute  values  of  the 
differences  between  the  scaled  frequencies  of  a given  case 
and  the  corresponding  experimental  frequencies  for  in-plane 
and  out-of-plane  modes  respectively.  The  calculated 
intensities  of  a given  method  were  multiplied  by  the 
corresponding  absorbance-intensity  conversion  constant  (as 
was  done  in  Section  3.4).  These  constants  ace  labelled  AICIP 
and  AICOF,  which  correspond  to  in-plane  and  out-of-plane 
nodes  respectively.  They  are  given  in  Table  5-21.  These 
converted  R.A.s  are  compared  to  the  experimental  R.A.s  of 
Table  3-9.  The  SUP  and  SIOF  values  are  the  mean  of  the 
absolute  values  of  the  deviations  of  the  calculated  R.A.s 
from  the  corresponding  experimental  R.A.s  for  the  in-plane 
and  out-of-plane  modes  respectively.  The  calculated  R.A.s 
for  4-31G,  SC,  HC,  and  STO-3G  cases  involved  the  use  of  the 
4-31G  APTs  (Table  3-5).  For  LC  case  the  LC  APTs  were  used. 
For  the  3GG,  XLl , and  XL2  the  3GG,  XC.1 , and  XL2  APTs 
rsspectively  were  used. 

Table  5-21  shows  that  the  best  approximation  to 
experimental  results  came  from  the  4-31G  case.  This  was 
followed  by  SC  and  HC  methods.  These  results  were  the 
expected  ones.  Some  randomness  is  observed  in  the  analysis 
of  the  calculations  of  the  out-of-plane  results.  These 
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15.1 
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577 

577 
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methods  were  all  better  than  the  ST0-3G  case.  The  results 
for  the  lower  convergence  studies  {LC,  3GG,  XLl,  and  XL2) 
did  parallel  to  a large  extent  the  observations  made  in  the 
previous  sections,  thus  giving  them  more  validity,  it  can  be 
said  that  to  some  extent  the  closer  the  predictions  of  a 
given  method  ace  to  the  4-31G  calculated  results  of  Chapter 
3,  the  closer  these  predictions  are  expected  to  be  to 
experimental  results.  It  Is  interesting  to  look  at  the 
absorbance-intensity  conversion  constants  of  Table  5-21.  The 
AICIP  values  are  rather  similar  for  all  the  methods  analysed 
here.  The  sane  is  true  for  the  AICOF  values.  The  average 
ratio  of  the  AICIP  to  Che  AICOP  value  is  1.5  for  a given 
method.  This  is  in  line  with  ocher  observations  where  it  was 
found  that  when  the  calculated  out-of-plane  intensities  were 
compared  with  experimental  results  the  calculated  values 
were  much  larger  for  out-of-plane  nodes  relative  to 


experiment 


1-plane  nodes. 


CHAPTER  6 

PREDICTION  OF  THE  IR  SPECTRUH  OF  2-PYRIMlDINONE 

This  chapter  is  concerned  with  the  prediction  of  the 
IR  spectruiQ  of  2-pyrinidinone . As  seen  in  Chapter  1 there 
was  no  experimental  spectrum  for  noninteracting,  isolated 
2-pyrinidinone  molecules.  In  the  gas  phase  and  in  dilute 
matrices  the  only  tautomeric  form  present  was  2-pyrinidinol . 
On  the  other  hand  the  only  tautomeric  form  present  in  the 
solid  phase  and  in  polar  solutions  was  2-pyrinidinone.  The 
IR  spectra  of  the  latter  system  corresponded  to  molecular 
aggregates  of  various  types.  This  fact  in  turn  creates  great 
difficulties  in  assigning  fundamental  bands.  Hence  the 
predicted  spectra  cannot  be  calibrated  by  comparison  with 
experimental  results. 

This  problem  for  2-pyclmidinone  is  alleviated  by  the 
fact  that  there  have  been  a great  many  experimental  and 
calculational  studies  of  the  infrared  spectra  of 
uracil . Based  on  the  similarity  in  structure 

between  the  diketo  form  of  uracil  {Figure  1-i)  and  that  of 
2-pytlmidinane  we  expect  to  find  similarity  between 

calculated  IR  spectra  for  these  two  molecules,  in  this 
chapter  we  first  examine  the  IR  spectra  of  uracil  to 
calibrate  our  calculations.  This  is  followed  by  the 
calculation  of  the  IR  spectrum  of  2-pyrimidinone . The 
chapter  is  divided  in  two  sections:  conclusions  from  uracil 
studies,  and  IR  spectrum  of  2-pyrimidinone. 

148 


6 . 1 Conclusions  from  Uracil  Studies 
t^ishiraura  et  calculated  Che  fundanental  in-plane 

vibcaCional  frequencies  of  Che  diketo  tautoner  of  uracil 
with  a "near  4-31C"  force  field.  This  force  field  was 
calculated  by  usinq  Equation  1-3,  where  was  the  force 

field  of  uracil  in  S coordinates  calculated  at  the  HF  5CF 
level  with  an  STO-3G  basis  set,  C was  a diagonal  matrix 
containing  4-31G/STO-3G  scaling  constants,  and  F was  the 
"near  4-31G”  force  field  in  S coordinates.  These  symmetrised 
internal  coordinates  are  given  in  Table  6-1.  The  conversion 
constants  4-31G/STO-3G  are  given  in  Table  6-2  under 
label.  This  F^^  force  field  was  calculated  using  the 
optimised  geometry  of  uracil  at  the  HF  SCF  level  with  the 
STO-3G  basis  set.  The  4-31G/5TO-3G  scaling  constants  were 
determined  from  calculations  of  force  constants  of  some 
smaller  model  compounds  similar  to  uracil  at  the  HF  SCF 
level  with  both  4-31G  and  STO-3G  basis  sets.  These  smaller 
molecules  were  "fragments"  of  uracil:  urea,  acrolein,  and 

formamide . 


repeat  the  calculation  of  the  in-plane  band  frequencies,  and 
PEDs,  and  extended  the  calculations  to  predict  also  the 
intensities  of  uracil.  The  APTs  needed  to  calculate  the 
intensities  were  determined  at  the  HF  SCF  level  with  a 4-31G 
basis  set  using  the  GAUSSIAN  76  program^^  at  a convergence 
of  5 * 10”^  and  at  an  integral  threshold  of  10"®.  The 
reference  state  used  in  the  calculation  of  these  apts  was 
the  geometry  reported  by  Nishimura  et  al.^®  based  on  x-ray 


structure,  atom  numbering  end  in-plane 
internal  coordinates  of  uracil  used 
in  this  work.  The  out-ol-plane  internal 
coordinates  are  the  out-of-plane  bends  and 
torsions  given  in  Table  6-1. 
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TABLE  6-1 

STHHETRIZED  INTERNAL  COORDINATES  TOR  URACIL 
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THE  4-31G/STO-3G 


SCALING  CONSTANTS  F 
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crysCallogcaphlc  data  of  1-cyclohexylucaci 1 [with  a few 
modifications  for  NH  and  CH  bonds].  Their  calculations  were 
in  essence  the  "5C  method*’  reported  in  Chapter  4.  They  used 
their  calculationa  to  assign  the  in-plane  experimental  IR 
fundamental  bands  of  uracil  with  considerable  success.  The 
experimental  IR  spectrum  was  obtained  from  investigations  of 
uracil  in  a dilute  Ar  matrix  at  10°  In  this  matrix 

uracil  was  found^*^  to  be  predominantly  in  Che  isolated 
diketo  tautomeric  form  (Figure  1-11.  These  studies  were 
reanalyxed  and  summarized  by  Scott-Lebron  in  his 

dissertation,^^ 

He  decided  to  reexamine  these  calculations  in 

accordance  with  our  procedures.  The  4-31G/STO-3G  scaling 
constants  and  APTs  of  uracil  wets  calculated  by  using  two 
sets  of  ^ initio  calculations.  One  set  was  calculated  at 
the  RHF  SCF  level  with  an  STO-3G  basis  set  using  subroutine 
FRBQ  of  GAUSSIAN  82^^  program  (Chapter  21,  at  the  geometry 
optimized  at  the  same  level  of  theory  (Nishimuca  et  al.^^) 
as  reference  state,  to  obtain  the  force  field.  He  refer  to 
this  force  field  as  the  "lew  field"  of  uracil.  The  other  set 
was  composed  of  the  diagonal  elements  of  the  force  field  of 
uracil  in  H coordinates  evaluated  at  HP  SCF  level  with  a 
4-31G  basis  set  using  the  GAUSIAN  78^^  program  (5  * 10  ^ 
convergence,  10~°  integral  threshold)  with  the  crystal 
geometry^^  as  reference  state.  He  refer  to  this  set  as  the 
"F  constants"  of  uracil.  Using  the  "low  field"  and  the  "F 
constants"  the  4-21G/5TO-3G  scaling  constants  were 
calculated  by  the  methods  described  in  the  discussion  of  the 
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wc  and  MCD  methoda  of  Chapter  4.  The  scaling  constants  for 
the  HC  work  are  not  reported  In  the  present  work  (they  can 
be  calculated  using  Equations  2-6  through  2-11).  The  scaling 
constants  tor  WCD  work  ace  reported  under  c^^^j  in  Table  6-2. 
The  APIs  of  uracil  were  calculated  at  HP  SCF  level  (with  a 
4-31G  basis  set)  using  GAUSSIAN  76®^  program  (at  5 • 10“® 
convergence!  10  integral  threshold)  using  the  crystal 
gaonetry^^  as  reference  state.  He  refer  to  them  as  the 
"crystal  APTs"  of  uracil.  The  in-plane  portion  of  the  "low 
field"  of  uracil  was  the  same  as  the  force  field  of 
Nishinura  et  al.^*^  The  out-of-plane  portion  of  the  "low 
field"  of  uracil  is  given  in  Table  6-3.  The  "crystal  APTs" 

Tables  6-4  through  6-6  show  the  calculated 

frequencies,  FEDs,  and  intensities  of  uracil.  The  ntrh 
results  in  these  tables  were  calculated  using  the 
scaling  constants  (Table  6-2),  the  "low  field,"  and  the 
"crystal  APTs"  of  uracil  by  the  "SC  method"  of  Chapter  4, 
i.e.,  the  same  way  as  the  uracil  calculations  of  Chin  et 
al.^^  The  (4C  results  were  calculated  using  the  4-31G/STO-3G 
scaling  constants  in  H coordinate  system  (not  reported  in 
the  present  work),  the  "low  field,"  and  the  "crystal  APTs" 
by  the  "(4C  method"  of  Chapter  4.  The  WCP  results  were 
calculated  using  the  4-31G/STO-3G  scaling  constants 
(Table  6-2),  the  "low  field,"  and  the  "crystal  APTs"  by  the 
"HCD  method'  of  Chapter  4.  Finally  the  STO-3G  calculations 
were  made  using  the  unsealed  "low  field"  and  the  "crystal 
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PREDICTED  BAND  FREQUENCIES  OF  URACIL  FROH 
THE  FORCE  FIELDS  OF  VARIOUS  HEIHODS 
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344  !6 


48.0 


; The  Creguencie 


methoas  are 


1S9 


PREDICTED  POTENTIAL  ENERGY  DISTRIBUTION 
FOR  URACIL  IR  FUNDAKENTAL  BANDS  BY  VARIOUS  METHODS 


1 100S7 

2 100S9 

3 91S11 

4 91S12 

5 S3S8 

6 S4S10 


S355 

loss 

16S10 


10S2 


99S7 

99S9 

43511 
SSS12 
S6S11 

43512 


30S10 

21SS 

39S10 


SPSS 


17S8 

IISIS 

2SSia 

isse 

Z0S16 

Z9S18 

I6S16 


11 


13517 

13518 
2481S 


30S21 

15S1 


10027 

10089 

22S11 

77S12 

78811 

22812 

S6S8 

11S18 

12S8 

S4810 


2286 

27S16 


10087 

10089 

90511 

90512 
72S8 

69S10 


S5S5 

12810 


22SS 


38816 


1981S 

38S21 


26S21  lesis 

13S17 


12 


21S3 

21S18 

19820 


14S16 

21S20 

43S21 

17S5 

3786 

11816 

24820 


11S2 

28S3 

16S18 

23820 

13816 

24820 

44821 


13SS 

37S6 

14816 

20820 


10S3 

25S16 

12S17 

13S2 

2783 

10818 

22S20 

12816 


148S 

4086 

21820 


11S3 

12816 

12817 

2483 

1084 

17818 


42S21 


12S5 


19820 

7781 


2SS4 


TABLE  6-5  continued 


47S1S 

14S19 

16S14 


24 


26 


28 


76S13 


11S3  26S4  1083 

2SS4  12513  2954 

12S13  10S13 

S7S15  SISIS  81S1S 


39S17  15S14  28814 

32S19  11S15  30S17 

33S17  20S19 

20S19 

70514  66S14  S0814 

16S15  14S1S  ISSIS 

12S17 


13S30 

67S22 

34S28 

15523 

41526 

15S27 

29S26 


14S23 

52S2S 

13S27 


68S24 

34S25 


21S22 

69S23 

13S38 

24S22 

16S26 


3SS19 

27S29 

91S30 


65S22 

33S28 

12523 


37S19 

36829 

S3S30 

20S28 

64S29 

21S30 

46S22 

46828 

91526 


18828 


24S26 

95527 

50S2S 

47826 


12S22 

52824 

37825 


17824 

39825 

13828 

50824 

59825 


100823  103823 


60522  55522 

S0S24  45524 


(Table  6-i: 
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Table  6-7  is  a sumnary  of  the  statistical  analysis  for 
the  pcedictions  for  these  different  methods.  The  CALC  set 
examines  errors  in  values  calculated  from  each  set  compared 
with  the  same  property  for  NTKM  predictions  of  Tables  6-4 
through  6-6.  and  values  of  Table  6-2.  The  £XPTL  set 

examines  differences  between  calculated  properties  and  the 
experimental  results.  The  experimental  results  are  the 
frequencies  and  relative  absorbances  of  the  experimental  IR 
spectrum  of  uracil  {from  studles^"^  in  a dilute  Ar  matrix  at 
low  temperatures}  as  analyzed  by  Soott-Lebron  These 
experimental  results  ace  plotted  schematically  in  Figures 
6-2  and  6-3  under  EXPTL  label.  In  order  to  measure  the 
statistical  parameters  of  the  EXPTL  set  of  Table  6-7,  the 
frequencies  predicted  by  a given  method  (Table  6-4)  were 
multiplied  by  0.9,  and  the  intensities  (Table  6-6)  were 
multiplied  by  0.603  (the  absorbance-intensity  conversion 
constant.  Figures  6-2  and  6-3  show  the  scaled  frequencies 
and  converted  intensities  (R.A.s)  calculated  by  the  various 
methods . 

The  CALC  set  of  Table  6-7  indicates  that  results  from 
the  "WCD  method"  were  closest  to  NTKH  pcedictions.  The  SXPTL 
set  of  this  table  indicates  that  NTKK  pcedictions  were 
closest  to  experimental  results,  followed  by  results  from 
the  "WCD  method."  Hence  it  could  be  said  that  the  closer  a 
method  was  in  reproducing  the  NTKK  predictions,  the  better 
the  prediction  when  compared  to  experimental  results.  From 
these  cough  trends  one  can  ma)ie  various  decisions.  The  most 
important  one  is  the  selection  of  NTKH  as  the  method 


which  imitaCes 


experimental  resulte.  Concerning 


Nishioiura  et  al. 
accurately  from 


ODie  interesting  conclusions  can  be  drawn, 
calculated  the  ^fj^KW  (Table  6-2) 

molecules  (fragments  of  uracil).  On 
values  (Table  6-2)  were  calculated 


from  HF  SCF  level  studies  with 
using  uracil  at  all  times.  Thu 
values  should  be  better  than 


4-llG  and  ST0-3G  basis  sets 
one  might  expect  that 


molecule.  As  seen  in  previous  chapters  this  would  probably 
be  true  if  the  "F  constants"  of  uracil  had  been  evaluated 
at  higher  convergence  (say  10~  instead  of  5 * 10  and 
with  the  optimised  geometry  (HF  SCF  with  a 4-31G  basis  set) 
instead  of  the  crystal  geometry^^  of  uracil.  The  case  is 
that  from  Table  6-7  it  can  be  concluded  Chat  were  the 


The  predicted  out-of-plane  frequencies.  PEDs,  and 
intensities  of  uracil  are  reported  in  Tables  6-4  through  6-6 
for  the  HC,  WCD,  and  STO-3G  cases.  They  are  not  given  for 
HTKK  since  Hishimura  et  al.^^  did  not  report  out-of-plane 
C^txm  values.  The  out-of-plane  predictions  from  these 
methods  are  not  statistically  analysed  in  this  study.  The 
major  reason  for  this  is  the  fact  that  there  were  many 
controversies  in  the  literature  about  the  assignments  of 
experimental  out-of-plane  modes  of  matrix  isolated 
uracil . Other  reasons  stem  from  the  statistical 
observations  on  out-of-plane  nodes  from  previous  chapters. 
In  Chapter  5 it  was  seen  that  in  some  situations  it  is  not 
advisable  to  scale  out-of-plane  force  fields  when 


rigure  6-2.  Conparison  between  the  calculated  IR  in-plane 
fundamental  inodes  ul  through  u6  fcon 
various  scaling  methods  and  the  corresponding 
fundamentals  from  experimental  results.  The 
experimental  results,  EXPT,  are  deduced  from 
Reference  61.  The  methods  of  calculation  are 
given  In  Chapter  6. 
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Wov«nuniD«ri  (c(n~'l 


3700 


fundamental  nodes  v7  through  v21  Cron 
various  scaling  methods  and  the  corresponding 
fundamentals  from  experimental  results.  The 

Reference  61.  The  methods  of  calculation  are 


1 > 

tA 

■4i 

T'%,’'  t’ 

" ?l  r ■K.i  1'^' , 

2000  IK» 

1000  sx 

16B 


STATISTICAL  ANALYSIS  OF  THE  PBRFOBHANCB  OF  VARIOUS  METHODS 
IN  IMITATING  THE  IN-PLANE  NTKM  AND  EXPERIMENTAL  RESULTS 


27.8  66.2 


32.7  38.6 


The  othec'^tecitis*ace  onitlesB?  The^terms^are  Hfined^ln^^’ 
Section  6.1. 


Uts. 
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convergence  ie  low.  See  Coc  inetence  the  LC  case  in  Table 
5-9.  These  situations  deteriorated  more  when  besides  low 
convergence,  the  reference  state  was  other  than  the 
optimized  geometry  (at  same  level  of  theory  as  the  one  used 
in  the  evaluation  of  force  field).  See  for  instance  the  STOP 
values  of  Tables  5-20  and  5-21.  Based  on  these  observations 
the  scaling  constants  evaluated  from  the  "F  constants"  are 
not  to  be  trusted  (low  convergence  and  crystal  geooietcy^^  as 
reference  state).  To  choose  the  best  set  of  out-cf-plane 

choice  being  the  STO-9G  case. 

6 . 2 IR  Spectrum  of  2-Pyrloiidinone 
The  complete  force  field  (in-plane  and  out-of-plane) 
for  2-pyrinidinone  was  evaluated  at  the  HF  SCF  level  with  an 
STO-3S  basis  set  using  the  FREQ  subroutine  of  GAUSSIAN  62^^ 
(Section  2. 1.1. 3).  The  reference  state  was  the  optimized 
geometry  at  this  sane  level  of  theory  (Section  2. 1.1.1).  The 
optimized  geometry  is  referred  to  as  3GG  of  2-pyrimidinone 
and  is  given  in  Table  6-6  (defined  in  Figure  6-4).  The  force 
field  is  referred  to  as  the  5TO-3G  force  field  of  this 
molecule  and  is  given  in  Table  6-9  in  S coordinates.  These 
symmetrized  internal  coordinates  are  given  in  Table  6-10. 

Ab  initio  calculations  of  diagonal  force  constants  and 
atomic  polar  tensors  of  2-pyrimidinone  were  carried  out  with 
GAUSSIAN  76^^  (at  5 * 10  ^ convergence,  10  ^ integral 

threshold)  at  two  reference  states.  These  reference  states 
were  the  previously  mentioned  3GO  and  another  labelled  as 
XL.  The  XL  geometry  was  obtained  from  x-ray  crystallographic 


Structuce,  atooi  nunbeclng  and  In-plana 
internal  coordinates  of  2-pyclnldinone  used 
in  this  work.  The  out-of-plane  internal 
coordinates  are  the  out-of-plane  bends  and 
torsions  given  in  Table  6-1. 


m 
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TKB  ST0-3G  0?T1HIZED  3GG  AND  CRYSTAL  STRUCTUBB  XL 
OF  2-FYRIHlDlNONE 


1.43S 


1.331 


1.219 

1.092 

1.077 


1.304 


1.355 

1.353 


1.244 

i.oao 

1.080 


12  123.1 


16  117.3 


19  123.0 

20  117.3 

21  122.4 

22  115.9 

23  119.9 

25  116.2 

26  120.3 


118.3 

118.1 


121.9 

121.6 

119.2 

116.4 

118.4 


113.7 

123.9 


TABLE  6-9 

FORCE  CONSTANTS  FOR  2-PYRIHlOINONE 
IN  SYHMETRIZSO  INTERNAL  COORDINATES 
EVALOATED  AT  THE  HF  SCF  LEVEL  WITH  AH  STO-3G  BASIS  SET 
g^a,b.e.d  S3  S4  S5  S6 


81 

S2 


7.23498 

.50644 

31723 


85 


-.41787 

.61717 

-.05468 

1.14269 

-.04685 

.02036 

-.06269 


21012 

05624 

02987 

51603 

00231 


S19 


-.04109 


6.11104 

.32673 

-.36551 

.37600 

-.30085 

-.00577 

.92016 

.02753 

-.02748 

.01233 

.28030 

-.33178 

-.00023 
. 39842 
-.03796 
.03122 
-.01373 


12770 

91375 

46001 

49190 

06699 

19571 

17022 

01933 

03274 


-.53746 

-.04014 

-.05208 

-.35736 

.01069 

-.00857 


6.89451 

.86423 

-.64967 

.05703 

-.05268 

-.02275 


37760 

62640 


-.16209 

-.01427 


02463 

03793 

06295 

27396 

06014 

07324 

15639 

00305 


06964 

01330 

37953 

04200 

00809 

06918 

00118 

22477 


6.70424 

.00243 

-.01995 

-.01627 

.02849 

-.19624 


33429 

19572 

15330 

02175 


TABLE  6-9 


S7  SB  S9  SIO  8U  S12 


SB 

S9 


512 

513 


BS04S  16.4B912 
00021  .01410 
0037$  -.01503 
000B2  .01060 
07967  .29278 
07387  .20326 
01904  -.33030 
01539  -.00664 
01B23  -.08925 
00110  -.01523 
00437  .00221 
02600  .02264 


7.36213 

-.00044 

-.00136 

.07641 

-.08181 

.02592 

-.00139 

-.00836 

-.01117 

.01734 

-.00559 


01693 

09091 

03906 


00176 

00873 

01507 

01505 

01457 


7.54319 

.06451  1.68422 

.01617  .03597 

.06567  .02032 

.02175  .00708 

-.01470  -.02402 

.00435  -.00926 

-.01599  .00429 

.00742  .02884 


S13  S14  515 


513  1.55059 

514  -.07313 


S17  .04017 

SIB  -.09364 


S19 


10019 


1.63200 

.13406 

-.02279 

.08272 

-.05930 

-.02133 


.66792 

.05793 

-.00561 

.01592 

.00640 


1.20466 

-.02565 

-.00757 

-.02758 


.77415 

.01349  .64601 

-.01691  .01104 


.69706 
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.04146  .1« 

.02910  -.03 


5 .75501 

2 -.05093  .66475 

5 -.00668  -.02830 


'“s  ■ 
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2-pyRIKIDlNONE 


(l/6)”'’(6B12-&ent6ei4'6SlS+aS16-AS171 
(1/12)® • ^i24B12-4613-i614*2fl61S-aei6-aB17) 
(1/2H  AB13-AS144AB16-ael7) 
(l/2)®-^(461B-aB23| 

(l/2)®-^lflS19-aS24. 

(l/2}®'®Ifl620-4825 
(l/2)®'®lfl821-flS26 
(l/2)®'^(4622-6827 


rNlC2 

rC2H3 

rN3C4 

rC4CS 

CC5C6 
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TABLE  6-10  continued 


Out-of-plang*' 


820  - (i/ei^'^iaTi-ira+ATi-aT^taTS-ate)  - ns 

521  - (i/2)(aTi-at3+aT4-a-tS}  - tiie 

522  - (i/i2)®'®(-aTi»2aT2-aT3-aT«+2aTS-aT«|  - xiiie 

823  - ArH7  - rH7 

524  - AvOS  - t08 

525  • BtH9  - tH9 

526  • arKlO  - yHIO 


®The  nomenclature  and  numberinq  according  to  Figure  6-4. 
The  recommendations  of  Fulay  et  al.  [ReCecence  11)  are 

^The  atoms  involved  in  torsions  ace:  Tle6123,  x2el234, 
X3-2345,  X4-34S6,  xS-4561,  and  T6-S612. 
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data^  (except  for  NH  and  CH  bonde  which  wece  modified  as 
decided  by  Chin  et  al.^®).  The  XL  georaetty  is  given  in  Table 
6-a  (accocding  to  Figuce  6-4).  The  3GS  foice  constants  and 
3GG  APIs  were  evaluated  with  3GG  as  reference  state.  The  XL 
force  constants  and  XL  APTs  were  evaluated  with  XL  as 
reference  state. 

By  the  "WCD  method"  of  Chapter  4 the  3GG  diagonal 
force  constants  and  the  STO-3G  force  field  of  2-pyrinidinone 
were  used  to  obtain  the  scaling  constants  Cjgj.  (Table  6-11) 
and  the  3GG  "near  4-31G"  force  field.  Similarly  Che  XL 
diagonal  force  constants  and  the  STO-3G  force  field  of 
2-pyrimidinone  were  used  to  obtain  Che  values  (Table 
6-11)  and  the  XL  "near  4-31G"  force  field.  The  3GG  and  XL 
APTs  are  given  in  Tables  €-12  and  6-13  respectively 
according  to  Che  orientations  of  Tables  6-14  and  6-15.  The 
3GG  frequencies,  PEDs,  and  intensities  of  Tables  6-16 
through  6-18  were  obtained  from  the  3GG  "near  4-31G"  force 
field  and  the  3GG  APTs.  The  XL  parameters  in  these  tables 
were  obtained  from  the  XL  "near  4-31G"  force  field  and  the 

Tables  6-16  through  6-lB  also  give  NTKH  and  STO-3G 
results.  The  NTEH  frequencies  and  normal  modes  were  obtained 
by  using  a "near  4-31G"  force  field  evaluated  by  the  "SC 
method"  of  Chapter  4.  This  "near  4-31G"  force  field  was 
evaluated  from  Che  scaling  constants  of  Nishimura  et  al.,®® 
^NTKM  Table  6-11,  and  the  STO-3G  force  field  of 

2-pyri»idinone  (Table  6-9).  The  NTKH  intensities  were 
obtained  by  using  this  "near  4-31G"  force  field  and  the  XL 


1-31G/ST0-; 


3G  SCALING  CONSTANTS  FOR  2-PTRlHIDINONE 
FROM  VARIOUS  METHODS 


10 

11 

n 


16 

17 

18 


82 

82 


6929 

6586 

7941 

7489 

8079 

7211 


7530 

7666 

7615 


7056 

7201 


6726 

7625 

7797 

7701 


.9645 

1.1364 

.8618 

.9611 

.9090 

[8470 


.8511 

.8409 

.8520 

.9645 

1.0153 

1.0169 


1.1025 

.0484 


8907 


.7653 

1.2021 


.8765 

1.0687 


5760 

3933 

7078 


!tho3‘ 


detaili 


-i-  ■#<  ,J.. 

„i.. 

■‘t”  -I-  -.J„, 

TiS!"’ 

-r  .J„ 

-F  ” ,..1. 

#'  -J,„ 

-r'i  ■#"  ,..L 

.70653  -.78716  0 

-.06643  .93014  0 

GsoneTRY 


TABLE  6-13 

ATOMIC  POLAR  TENSORS 
FOR  2-PVBiniDINONE  HITS  CRYSTAL  XL 


-.081S4  .905S6 


^APTs  ace  detecoined  at  HF  5CF  level  with  4-31G  basis  set 
using  crystal  geonetcy  (Table  6-10)  as  ceCecence  state. 


systera  of  orientation  is  the  principal  cartesian  coordinates 
of  Table  3-6. 


X® 


-.200061 

-1.000213 

-.271430 

1.02149S 

1.033102 

1.192479 

-.696748 

-2.217638 

1.S49207 

2.910150 

1.680693 


3.21680641  ( 10‘*  CG,  lO*'^ 
1.04140900  (5  * 10'^  CG, 


DIPOLE  (PEBYES)*’ 


• 7.1117  (5  * 10“^  CG,  10'®  IT,  4-31G  BS)^ 


PAL  CARTESIAN  COORDINATES.  GROUND  STATE  ENERGY, 
AND  DIPOLE  nOHENT  OF  2-PYBIHlDINONE 
AT  THE  XL  GEOKETRY 


C5 

C6 

HI 


1.131739 

-.030594 

-1.232909 

-1.195651 

-.025936 

1.152426 

2.025603 

.026201 


-.207724  0 

-.949650  0 

-.201013  0 

1.03B705  0 

K144469  0 

-.697512  0 

-2.192353  0 


HIO  -.018713  2.893391  0 

Hll  2.123810  1.616500  0 


ENERGY  lAU)^ 

- -337.04532768  IS  • 10“^ 


DIPOLE  (DEBYES)** 

- 7.8810  (5  * 10"®  CO.  10*®,  4-31G  BS) 


^Principal  cartesian  coordinates  are  in  A,  The  crystal  XLl 
geometry  is  shovn  in  Table  6-S.  The  atom  numbering  is 
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APTft  (Table  6-12).  The  ST0-3G  frequencies,  PEDe,  and 
Intensities  were  obtained  Ciom  the  unseated  sto-3G  force 
field  of  2-pyribidinone  (Table  6-9)  and  the  XL  APTs  (Table 

Table  6-19  sumnarizes  the  statistical  analysis  for  the 
In-plane  predictions  of  the  various  methods  of  Tables  6-16 
through  6-18.  The  errors  are  measured  with  respect  to  NTXn 
predictions  for  2-pyrinidinone . The  MOIP  and  HIP  values  are 
measured  with  respect  to  values  o£  Table  6-11.  In  this 
Table  one  can  see  that  XL  case  imitated  best  the  NTKH 
results,  followed  by  6TO-3C  and  3GG  cases.  Although  we  do 
not  have  an  experimental  spectrum  of  isolated  2-pyrimldinone 
we  believe  from  the  uracil  observations  that  NTRH  is  the 
closest  to  experimental  results  followed  by  XL.  Thus  we 
choose  the  ntkm  predictions  of  Tables  6-16  through  6-18  as 
best  for  the  in-plane  fundamental  modes  of  the  experimental 
spectrum  of  isolated  2-pyrimidinone . As  in  the  observations 
on  uracil  one  is  confronted  with  a choice  of  scaling 
constants  (Table  6-11)  from  which  values  are  chosen  as 
best  followed  by  in-plane  constants. 

Arguing  from  nonstastlsti cal  angles  one  can  arrive  at 
similar  conclusions.  The  "near  4-31G"  force  fields  of  ntkk 
cases  for  uracil  and  2-pyrimidinone  cases  were  obtained  by 
scaling  the  force  field  (at  the  same  level  of  theory)  of  two 
closely  related  molecules  by  the  same  set  of  conversion 
constants  APTs  for  both  molecules  were 
evaluated  at  the  same  level  of  theory  with  similar  crystal 
structures  as  reference  states.  Thus  if  the  ntkh  case  for 
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P8S0ICTED  POTENTIAL  ENERSY  DISTRIBUTION  FOR  2-FYRIHlDINOL 
IR  FUNDAMENTAL  BANDS  FROM  FORCE  FIELDS  OF  VARIOUS  METHODS 


1 100S7 

2 92S10 

3 92S11 


66SS 

23S3 

3SS5 

10S19 

2653 

20SS 

1253 
30S6 

16515 

16S18 

11S19 

1254 
29S15 

1SS19 

10S2 

21S3 

47S17 

3181 

ISSIS 

12S17 

1052 

2SS16 

40S19 


15  S2S12 

16  2881 
S0S2 


100S7  100S7  100S7 

94510  79S1D  92510 

19S11 

94511  19510  92511 

9959  98S9  9959 

8158  7158  7958 

1853  2053  2953 

S4SS  46S5  39SS 


50S3 

1154 

2556 

27S7 

11519 


3983 

1254 


3156 

25515 


4083 

2755 


3156 

24S1S 
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PREDICTED  BAND  INTENSITIES  FOR  2-PTRiniDINONE 
FROK  VARIOUS  FORCE  FIELDS  AND  ATONIC  POLAR  TENSORS 
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ucacil  gave  successful  coc relations  with  in-plane 
experiaental  findings  so  should  the  NTKil  case  of 
2-pyriinidinone . 

Thus  one  can  be  confident  that  the  in-plane  HTKH 
predictions  for  2-pyriraidlnone  ace  quite  close  to  the 
expecinental  results.  As  far  as  the  out-of-plane  predictions 
are  concerned  we  do  not  have  a satisfactory  answer.  Based  on 
observations  of  the  out-of-plane  inodes  of  uracil,  we 
conclude  that  it  is  safer  to  use  the  STO-3G  out-of-plane 
predictions  of  Tables  6-16  through  6-18  as  the  best  set  of 
results.  The  out-of-plane  predictions  for  3GG  and  XL  cases 
are  reported  in  order  to  provide  researchers  with  more  tools 
at  their  disposal  but  not  much  more  can  be  said  about  them. 

we  conclude  this  chapter  by  presenting  Figure  6-5 
which  compares  Che  unsealed  predicted  spectra  of  the  two 
tautomers  of  2-oxopyr imidine . In  this  figure  Che  spectrum  of 

2- pyrlmidinol  is  obtained  from  the  unsealed  results  of  Table 

3- 7.  This  corresponds  to  our  beet  predicted  spectrum  for 
this  tautomer  at  Che  "full  4-31G"  study  of  Chapter  3.  The 
spectrum  of  2-pytimidinone  in  this  figure  depends  on  two 
sets  of  predictions.  The  in-plane  predictions  ace  from  the 
NTKH  results  and  Che  out-of-plane  from  Che  STG-3G  results 
(Tables  6-16  and  6-16).  This  is  our  predicted  best  spectrum 
for  2-pyrimidinone.  The  quality  of  the  latter  should  be 
inferior  to  the  2-pyrimidinol  spectrum. 
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STATISTICAL  ANALYSIS 
OF  THE  ABILITY  OF  THE  VARIOUS  METHODS 
IN  IMITATING  THE  NTKH  RESULTS  FOR  2-PYBIMIDINONE 
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Figure  6-5.  Simultaneous  display  of  best  calculations 

of  the  present  work  for  the  tR  spectra  of  the 
tautomers  of  2-oxopy rimldine . The  results 
for  2-pyrinidlnol  are  the  unsealed  calculations 
of  the  full  4-310  study  (Table  3-7).  The 
results  for  2-pyrimldlnone  correspond  to 
In-plane  NTRM  and  out-of-plane  ST0-3G  methods 
(Tables  6-16  and  6-16) . 
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CHAPTER  7 

SOHHARY  AND  COMMENTS 

The  calculations  of  the  IR  spectra  of  l-oxopyrinidine 
tautomers  are  quite  useful  for  studies  of  molecules  of 
similar  structure.  For  Instance  they  may  be  helpful  In 
analyeing  experimental  matrix  IP  spectra  of  nucleic  acid 
bases  which  exhibit  various  tautomeric  forms,  and  of  the 
halogenated  l-oxopytlmldlnes.®  These  results  can  be  further 
implemented  in  more  ambitious  aims  trying  to  Interpret  the 
spectra  of  nucleotides. 

The  calculated  spectra  are  helpful  in  order  to 
understand  better  the  experimental  results.  This  is 
particularly  true  for  large  molecules  with  low  symmetry  or 
short  lived  species  such  as  radicals  and  ions.  The  methods 
and  results  of  chapters  S and  S can  he  of  utility  In  further 
attempts  to  Improve  the  quality  of  calculated  spectra  by 
scaling  to  fit  other  calculations  or  experimental  results. 
The  mathematical  and  statistical  level  of  these  analyses 
clearly  needs  to  be  developed  more.  The  various  tables  in 
these  chapters  include  results  which  can  be  further  analysed 
in  the  future  when  more  sophisticated  approaches  are 
conceived . 

The  assignment  of  spectra  by  referring  to  calculated 
frequencies  only  is  demonstrated  to  be  insufficient.  This 
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approach  is  analogous  to  a mathematician  trying  to  £ind  the 
values  of  variables  from  n simultaneous  equations.  The 
predicted  normal  modes  and  the  predicted  band  intensities 
are  expected  to  be  helpful  in  these  correlations.  Person  and 
huBulat  are  developing  a method  analogous  to  the  FED  case 
whereby  the  intensities  are  described  in  terms  of  internal 
coordinates.^*^  Correlations  of  spectra  based  on  predicted 
frequencies  and  intensities  (ideally  simultaneous  instead  of 
independently) , in  conjunction  with  their  corresponding 
interpretations,  would  provide  the  mathematician  of  the 
analogy  with  note  relationships  that  would  help  him  to  solve 
his  simultaneous  equations  with  less  assumptions. 

Correlations  of  spectra  from  calculations  based  on  the 
usage  of  simple  basis  sets  such  as  4-31G  and  STO-30,  is 
recommended  as  a reasonable  system  for  Che  development  of 
more  sophisticated  mathematical  and  statistical  methods  of 
analysis.  This  Is  partly  based  on  the  idea  that  the 
wavefunctlon  of  the  molecule  is  better  understood  from  these 
basis  sets  than  from  experimental  facts. 

Various  observation  are  made  with  regards  to 
convergence,  integral  threshold,  and  reference  state.  It  is 
found  that  a desirable  level  of  precision  must  be  observed 
in  ab  initio  calculations.  The  choice  of  geometry  to  be  used 
at  different  times  in  the  calculations  is  a difficult 
problem.  This  "geometry  parameter"  enters  in  the 
calculations  at  various  points;  reference  states  for 
evaluation  of  force  constants,  APTs,  frequencies,  and 
intensities  of  IR  spectra.  Our  observations  indicate  that  a 


fore*  field  calculated  at  a given  level  of  theory  ahould  use 
for  a reference  state  the  geometry  optimized  at  the  sane 
level  of  theory.  Scaling  nethods,  including  other  nethods 
besides  the  presented  ones,  deserve  further  study. 

Due  to  the  high  cost  and  further  complications 
Involved  in  ab  Initio  calculations  of  IB  spectra, 
particularly  of  large  molecules,  various  researchers  are 
compiling  data  from  various  research  centers  as  well  as 
synchronizing  the  nature  of  these  calculations.  Along  these 
lines  we  have  calculated  4-31G/STO-3G  conversion  constants 
for  the  in-plane  nodes  of  pyridine,  using  our  results 
calculated  with  an  STO-3G  basis  set  together  with  the 
calculations  of  Wiberg  et  al.^^  who  used  a 4-31G  basis  set. 
These  conversion  constants  were  surprisingly  similar  to  the 
C constants  (Table  4-3)  for  2-pyrimidinoI  implying  some 
sense  of  universality  in  these  calculations.  Similar  work 
including  experimental  results  can  be  seen  in  the  work  of 
Pongor  et  al.^^  who  combined  data  on  the  out-of-plane 
nodes  of  py ridi  ne , or  in  the  conclusions  of  Harsanyi 

et  al."^^  who  combined  data^^'^^'^^  on  the  in-plane  and 
out-of-plane  nodes  of  uracil.  More  thorough  reviews  can  be 
seen  on  previously  indicated  references . It  is  in 
this  spirit  of  cooperation  that  most  of  the  work  of  Chapters 
4 and  5 was  carried  out. 

In  the  present  work  the  usage  of  ^ initio 
calculations  at  the  HP  SCF  level  with  basis  sets  such  as 
ST0-3G  and  4-31G  shows  to  be  helpful  in  the  Interpretation 
of  IB  spectra  (particularly  in  the  case  of  large  molecules 


wibh  low  eyminetry),  The  level  of  correlation  between  the 
presented  expe rinental  and  theoretical  results  is  good.  The 
objective  of  reliable  predictions  of  the  spectra  of  the 
tautoners  of  2-oxopyrimidine , and  of  the  interpretation  and 
correlation  with  experinental  matrix  observations,  is 
achieved  at  a reasonable  level  of  confidence.  The  more 
generalized  goals  of  investigating  the  quality  of  ^ initio 
calculations  of  IR  spectra  do  lead  to  enlightening  but 
certainly  not  conclusive  facts.  It  is  emphasized  that  these 
calculations  need  to  be  studied  at  greater  depth  in  the 
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